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Abstract 

The  feasibility  of  using  the  mass  ratio  Ndl46/Ndl45 
and  the  activity  ratio  Csl3VCsl37  as  off-site  Indicators 
of  bomb-grade  plutonium  production  in  commercial  light 
water  reactors  is  examined.  The  theoretical' basis  of  these 
ratios  as  on-site  indicators  of  fluence  is  developed 
and  experimental  validation  of  their  utility  as  such  is 
examined.  The  escape  rate  of  cesium  and  neodymium  is 
approximated  and  detection  limits  determined.  Although 
both  ratios  are  experimentally  validated  as  indicators 
of  fluence,  neither  cesium  nor  neodymium  escapes  a reactor 
in  great  enough  concentration  to  be  accurately  measured 
9 off-site . 
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ISOTOPIC  CORRELATION  TECHNIQUES 
AS  AN  OFF-SITE  REACTOR  MONITOR 

I.  Introduction 


The  United  States  would  like  to  export  nuclear 
reactor  technology  and  fuel  without  contributing  to 
the  proliferation  of  nuclear  weapons.  Prospective  third 


r. 


world  customers  do  not  possess  the  technology  to  convert 
commercial  uranium  fuel  to  bomb-grade  uranium;  this 
requires  an  inversion  of  the  ratio  of  U238/U235>  which 
is  97/3  for  a representative  light  water  reactor  (LWR) 
fuel.  However,  absorption  of  fast  neutrons  from  fission- 
ing U235  converts  U238  into  Pu239-  Pu239  has  nuclear 

properties  similar  to  those  of  U235;  it  can  be  used  as 
reactor  fuel  or  weapon  fuel.  The  continuous  conversion 
of  U238  to  Pu239  in  a power  reactor  burning  U235  is  the 
proliferation  concern. 

If  the  Pu239  Is  not  removed  from  the  reactor  core 
soon  after  its  formation,  it  too  absorbs  neutrons  and  Is 
converted  to  ?u2^0.  Substantial  amounts  of  Pu240  in  a 
mixture  of  the  two  plutonium  Isotopes  makes  the  plutonium 
unusable  in  a weapon;  moreover,  the  separation  of  these 


isotopes  is  even  more  difficult  than  the  enrichment  of 
uranium.  Economical  generation  of  power  requires  a long, 
sustained  fission  cycle;  production  of  bomb-grade  plutonium. 
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shutdown  for  Pu239  removal.  Obviously,  on-slte  inspection 
to  determine  length  of  fission  cycle  is  the  most  straight- 
forward method  for  determining  whether  a reactor  is  being 
used  for  production  of  power  or  Pu239  isolation. 

Another  method  for  determining  Pu239  isolation  in- 
volves examination  of  irradiated  fuel  elements.  Weitkamp 
has  demonstrated  a linear  relationship  between  U235 
depletion  and  the  ratio  Pu2^0/Pu239  (Ref  1:212).  Data 
generated  by  Origen,  the  Oak  Ridge  National  Laboratories 
Isotope  Generation  and  Depletion  Code,  verifies  this 
relationship;  U235  depletion  and  the  ratio  Pu2^0/Pu239 
have  a linear  correlation  factor  of  0.9995  (Ref  2:106). 
Determination  of  U235  depletion,  or  burnup,  of  a fuel 
element  indicates  the  relative  amounts  of  Pu240  and  Pu239; 
low  burnup  infers  Pu239  isolation.  One  method  of  deter- 
mining fuel  burnup  is  called  Isotopic  correlation  tech- 
niques (ICT);  It  is  based  on  the  fact  that  the  ratio  of 
certain  fission  product  nuclides  varies  systematically 
with  fuel  burnup.  Use  of  a ratio  eliminates  the  need 
for  knowledge  of  the  original  fuel  mass. 

ICT  have  demonstrated  their  utility  as  an  on-site 
safeguards  technique  for  determining  Pu239  isolation; 
they  are  also  used  at  fuel  recycling  facilities  to  verify 
fissile  content  of  incoming  shipments.  The  problem 
addressed  herein  is  the  feasibility  of  using  Isotopic 
correlation  techniques  to  monitor,  off-3ite,  for  production 
of  bomb-grade  plutonium  in  a power  reactor. 
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Goodwin,  dealing  with  the  same  prob  lem,  developed 
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a set  of  criteria  which  a nuclide  ratio  must  meet  if 
off-site  application  of  ICT  is  to  be  possible  (Ref  3:10). 

The  essential  elements  of  those  criteria  are: 

(1)  The  nuclides  must  be  isotopes  of  the  same 
element . 

(2)  This  element’s  chemical  properties  must  enable 
it  to  escape  from  the  fuel  and  reactor. 

(3)  The  istopic  ratio  must  be  measurable  off-site. 
Goodwin  predicts  that  in  1982,  79%  of  the  world’s  reactors 
will  be  LWR’s,  and  that  92%  of  the  world's  reactors  will 
be  using  slightly  enriched  (1-3%)  uranium  dioxide  (Ref 
3:13,15).  Only  LWR's  burning  slightly  enriched  uranium 
dioxide  will  be  considered. 

Isotopic  correlations  have  been  the  subject  of  both 
theoretical  and  experimental  Investigation  (Ref  4:426). 
Theoretical  Investigation  is  generally  carried  out  by 
performing  ac.curate  burnup  calculations,  from  which 
correlations  are  derived.  Experimental  investigation  is 
based  on  detailed  measurement  of  irradiated  fuel  isotopic 
composition,  from  which  correlations  are  derived. 

Analysis  of  the  off-site  application  of  ICT  began 
with  a literature  survey  of  theoretical  correlation  in- 
vestigations. Isotopic  ratios  thereby  identified  were 
evaluated  against  Goodwin's  criteria;  the  ratios  Csl3VCsl37 
and  Mdl46/Ndl45  were  Judged  suitable.  A literature  survey 
of  experimental  investigations  of  these  two  ratios  was 
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conducted  to  validate  them  as  indicators  of  fuel  depletion. 

( 

The  literature  revealed  a power  history  dependence  of 
the  cesium  ratio,  but  the  degree  of  this  dependence 
and  its  impact  on  the  utility  of  the  cesium  ratio  as  a 
burnup  indicator,  independent  of  knowledge  about  power 
history,  were  not  clear.  Calculations  were  performed  to 
determine  the  degree  and  impact  of  this  dependence.  The 
escape  mechanisms  from  fuel  to  environment  were  examined 
and  escape  rates  estimated.  These  escape  rates  were 
coupled  to  the  dilution  processes  operating  on  the  two 
reactor  effluents  considered,  stack  gas  and  coolant 
water,  and  limits  of  ratio  measurability  predicted. 

( 


II . Theoretical  Analysis  of  ICT 


The  evolution  of  isotopic  correlation  techniques 
began  with  W.  J.  Maeck's  proposal,  in  1965,  that  ratios 
of  fission  product  nuclides  be  used  for  determination 
of  nuclear  fuel  burnup  (Ref  5).  At  that  time,  determina- 
tion of  burnup  by  fission  product  analysis  was  based  on 
the  quantitative  measurement  of  a selected  fission  product 
nuclide.  Use  of  a ratio  instead  of  an  absolute  atom 
abundance  circumvented  the  requirements  for  complete 
sample  dissolution  and  volume  measurements. 

Maeck's  method  was  based  on  the  premise  that  for 
stable  fission  product  A,  with  large  capture  cross 
section,  and  stable  fission  product  A+l,  with  small 
capture  cross  section,  the  ratio  A+l/A  increases  ex- 
ponentially with  irradiation  time.  His  calculations 
indicated  that  the  ratios  Ndl44/Ndl43,  Kr84/Kr83,  and 
Xel32/Xel3I  were  suitable  for  burnup  determination. 

The  first  interest  in  radionuclide  ratios  as 
indicators  of  fuel  burnup  was  demonstrated  in  1968  by 
Hick  and  Lammer  (Ref  6).  They  developed  a computer  code, 
IRREL,  for  calculation  of  fission  product  activities  as 
a function  of  fluence  and  determined  that  the  activity 
ratio  Csl34/Csl37  was  an  effective  measure  of  fluence. 
Fluence  is  a measure  of  fuel  exposure;  fuel  depletion 
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Is  easily  calculated  from  the  relationships: 


ML  = 


255 


Mo?)  e^f) 


Uz35'  htPcz  r/oK/  — (4  - <y  ^ 


(i) 


(2) 


where 


A/co)  = initial  amount  of  U235 


2JS" 


A/  = amount  of  G235  at  fluence  (f) 

' 2ss 


0^  = absorption  cross  section  or  U235 

= fluence 

The  most  complete  compilation  of  isotopic  ratios  for 
LWR  fuels  is  that  of  the  Joint  Research  Center  of 
Euratom  (Ref  *4:427).  They  identify  the  following  ratios 
as  indicators  of  U235  depletion:  Csl34/Csl37,  Kr84/Kr86, 
Kr86/Kr83,  Kr84/Kr83,  Xel32/Xel31,  Xel34/Xel31,  Xel32/Xel34, 
Ndl46/Ndl45,  Ndl46/Ndl48,  Ndl48/Ndl45.  Two  ratios  of 
neodymium  isotopes  identified  by  Maeck  (Ref  5:15)  are 
ignored  in  the  preceding  compilation;  they  are  Ndl44/Ndl43 
and  Ndl46/Ndl43.  These  twelve  ratios  meet  the  first  of 
Goodwin's  criteria.  Their  escape  potential  and  off-site 
measurability  must  be  determined.  Since  escape  potential 
is  determined  by  chemical  characteristics,  use  of  same- 
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element  ratios  simplifies  escape-path  analysis.  The  im- 
portant determinants  of  escape  potential  are  state  (solid, 
liquid,  gas)  and  chemical  reactivity.  Escape  potential 
increases  as  state  changes  from  solid  to  gas  and  decreases 
as  reactivity  increases. 

Krypton  and  xenon  are  inert  gases;  their  escape 
potential  is  extremely  high.  However,  all  their  isotopes 
of  interest  herein  are  naturally  occurring.  Kr83  and 
Xel31  have  isomeric  states,  but  the  half-life  of  Kr83m  is 

_3 

less  than  two  hours  and  only  8E-3(i.e.  8 x 10  ) of  the  Xel31 

passes  through  the  isomeric  state  (Ref  7).  Disregarding 
these  isometric  states,  calculations  were  performed  to 
compare  expected  effluent  concentrations  of  the  three 
krypton  isotopes  and  the  three  xenon  isotopes  to  environ- 
mental background  concentrations.  Using  data  from  the  BWR 
at  Oyster  Creek,  N.  J.,  it  was  determined  that  the  stack 
gas  contained  concentrations  on  the  order  of  1E6  per  cc  of 
the  Krypton  isotopes  and  1E7  per  cc  of  the  Xenon  isotopes. 

The  atmospheric  background  concentrations  were  calculated 
to  be  on  the  order  of  1E12  per  cc  for  krypton  isotopes 
and  1E11  for  xenon  Isotopes.  These  calculations  are  con- 
tained in  Appendix  A.  Krypton  and  xenon  are  not  present 
iri  measurable  amount  in  coolant  water  effluent.  Clearly, 
krypton  and  xenon  do  not  meet  the  off-site  measurabili- 
ty criteria;  they  will  not  be  considered  further. 
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The  escape  potential  of  the  neodymium  Isotopes  would 
appear  to  be  low;  neodymium  is  a non-volatile  reactive 
rare-earth  metal.  Its  melting;  point  is  1024  degrees 
(centigrade)  and  its  boiling  point  is  3027  degrees  (Ref  8:B 
121).  Escape  data  kept  by  reactor  operators  and  agencies 
regulating  them  deal  only  with  radionuclides  (Ref  9).  Of 
the  neodymium  isotopes  of  concern  herein,  only  Ndl44  is 
radioactive,  but  its  half-life  of  5E15  years  causes  its 
radioactivity  to  be  commonly  disregarded.  As  demonstrated 
in  Appendix  A,  measured  concentrations  of  one  isotope 
can  be  used  to  infer  concentrations  of  other  isotopes  of 
that  element.  The  only  indication  of  neodymium  escape  found 
in  the  literature  was  a Ndl47  concentration  of  3E-7  Curies 
per  liter  of  coolant  in  a BWR  (Ref  10:8).  This  activity 
concentration  equates  to  a number  density  of  1.54E10  atoms 
of  Ndl47  per  liter.  Thermal-neutron  fission  of  U235  results 
in  the  following  yields: 

6.03%  Ndl43,  5.62%  Ndl44,  3-98%  Ndl45, 

3.07%  Ndl46 , 2.70%  Ndl47,  1.71%  Ndl48 
(Ref  11:98).  All  these  isotopes  could  be  expected  to  have 
been  present  in  the  coolant  with  concentrations  of  the  same 
order  of  magnitude  as  Ndl47. 

Although  1E10  atoms  per  liter  is  certainly  a measur- 
able amount,  neodymium's  ability  to  escape  from  the  coolant 
to  the  environment  is  uncertain.  If  neodymium  does  escape 
to  the  environment,  its  off-site  measurability  will  be  com- 
plicated by  the  natural  occurrence  of  all  Isotopes  of 
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interest  (Ref  8:B51).  Further  consideration  of  neodymium's 
escape  potential  and  off-site  measurability  is  deferred 
until  the  analysis  of  escape  mechanisms  is  completed.  Since 
tne  neodymium  ratios  are,  at  this  point,  of  questionable 
utility,  experimental  validation  for  only  one  of  the  ratios 
will  be  considered.  In  addition  to  being  cited  as  an 
indicator  of  U235  depletion,  the  ratio  Ndl46/Ndl^5  shows 
a direct  relationship  to  Pu2*IO  content  (Ref  1:212);  only 
this  neodymium  ratio  will  be  considered  further. 

The  chemical  characteristics  of  cesium  do  not  immedi- 
ately categorize  its  escape  potential;  although  volatile 
(boiling  point  of  690  degrees),  it  is  the  most  electroposi- 
tive (reactive)  element  (Ref  8:0106).  However,  one  environ- 
mental study  has  shown  that  0813*1  and  Csl37  can  be  found  in 
measurable  concentrations  in  both  the  stack  gas  and  discharge 
canal  of  a PWR  (Ref  12).  Since  neither  Csl3**  nor  Csl37 
occur  naturally  (Ref  8:B46),  off-site  measurability  will 
not  be  complicated  by  natural  background  concentrations. 

The  cesium  ratio  appears  to  be  the  best  candidate  for  off-site 
application. 

The  literature  survey  of  theoretical  investigations 
of  ICT  upon  which  the  preceding  selection  of  ratios  was 
based  revealed  three  other  items  of  interest.  The  first  of 
these  items  is  the  creation  by  the  International  Atomic 
Energy  Agency  (IAEA)  of  an  ICT  data  bank;  data  has  been 
collected  from  reprocessing  plants  and  analytical  labora- 
tories and  stored  with  a common  format  for  ease  of  subse- 
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quent  retrieval  (Ref  13).  The  ratio  chosen  as  an  on-site 
indicator  of  U235  depletion  was  Csl3^/Csl37-  A code  (ISOCORR) 
has  been  developed  which  performs  the  theoretical  analysis 
of  correlations  and  compares  the  calculated  data  with  the 
experimental  data  stored  in  this  bank(Ref  1*1). 

The  second  item  of  interest  is  a wide  concern  over 
the  sensitivity  of  theoretical  isotopic  correlations  to 
fission  product  nuclear  data  (FPND)  accuracies  (Ref  1,  15-17). 
Weitkamp  identifies  the  yield  of  Csl37  as  the  only  signifi- 
cant FPND  involved  in  the  production  scheme  of  Csl37  (Ref 
17:201).  He  apparently  disregards  the  short-lived  precursors 
of  Csl37  (1137  and  Xel37).  Weitkamp  lists  the  yield  and 
capture  cross  section  of  Csl33  and  the  capture  cross  section 
of  Csl3^  as  the  significant  FPND  in  the  production  scheme 
of  Csl3^;  again,  he  has  disregarded  the  precursors  of  Csl33. 

As  the  uncertainties  in  cross  sections  are  greater  than 
the  uncertainties  in  yields,  the  accuracy  of  the  cross 
sections  of  Csl33  and  Csl3^  is  the  limiting  factor  in  cal- 
culating a Csl3^/Csl37  ratio.  A 1 1 error  in  the  cross 
section  of  Csl33  causes  a 1%  error  in  the  amount  of  Csl3^ 

(Ref  17:200).  The  Important  FPND  in  the  calculation  of  the 
Ndl46/Ndl^5  ratio  are  capture  cross  sections  of  NdlM,  Ndl^5 
and  yields  of  Ndl^5,  Ndl46  (Ref  17:205). 

The  final  item  is  an  indirect  reference  to  the  power 
dependency  of  the  Csl3VCsl37  ratio  as  a measure  of  fluence; 
this  dependency  was  to  become  a concern  during  examination 
of  the  experimental  investigations  of  ICT.  Fig.  1 shows  the 


theoretical  variation  of  the  Csl3VCsl37  activity  ratio 
with  irradiation  time  and  neutron  flux  for  fission  of  U235. 
The  intersections  of  the  flux  and  time  lines  are  the 
plotted  points.  The  two  circled  points  both  represent 
the  same  fluence;  varying  the  flux  by  a factor  of  two 
changes  the  Csl3VCsl37  ratio  by  21%. 
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Pig.  1.  Csl3VCsl37  Activity  Ratio  at  Shutdown  Versus  Integrated 
Neutron  Flux  for  Different  values  of  Neutron  Flux  and  Irradiation 
Time  (Ref.  15:241) 


III.  Experimental  Validation  of  Cesium  and 
Neodymium  Ratios 

Selection  of  the  ratios  Csl3VCsl37  and  Ndl46/Ndl^5 
as  prospective  off-slte  Indicators  of  fuel  turnup  is  based 
on  their  theoretical  utility  as  on-site  turnup  Indicators; 
experimental  validation  of  their  on-site  utility  was  the 
subject  of  a second  literature  survey.  The  work  of  Hick 
and  Lammer  (Ref  18)  demonstrates  an  interesting  application 
of  both  the  theoretical  and  experimental  aspects  of  isotopic 
correlations;  they  combine  gamma-spectrometric  measure- 
ments and  theoretical  calculations  to  establish  a fission 
product  inventory  which  is  more  accurate  and  complete  than 
is  possible  with  either  experiment  or  calculation  alone. 
Their  code  (IRREL)  for  calculation  of  fission  product 
Inventory  requires  neutron  flux  Input  data  (value  and 
spectral  shape)  of  unavailable  accuracy.  On  the  other 
hand,  gamma-spectrometric  measurements  are  sufficiently 
accurate  for  only  a limited  set  of  fission  products;  due 
to  mutual  line  interference  or  absence  of  gamma  lines, 
many  isotopes  cannot  be  measured  accurately.  Hick  and 
Lammer  use  experimental  results  for  Cp13^  and  Csl.37  and 
a knowledge  of  the  irradiation  history  to  derive  proper 
neutron  flux  input  data  for  IRREL;  IRREL  then  is  able  to 
compute  a fission  product  inventory  of  any  desired  degree 
of  completeness. 

1 3 


The  relationship  between  fluence  and  fuel  depletion 


is  defined  in  the  preceding  chapter;  another  measure  of 
fuel  depletion,  or  burnup,  is  preferred  in  experimental 
work.  This  unit  is  the  megawatt-day  per  tonne  (MWD/T).  As 
the  name  suggests,  fuel  in  a reactor  operating  at  a specific 
power  level  of  one  megawatt  (thermal)  per  tonne  (1000  Kilo- 
grams) of  fuel  initially  present  will  be  depleted  by  one 
MWD/T  in  one  day.  The  fuel  is  considered  to  consist  of 
heavy  elements  only;  in  fuel  elements  of  uranium  dioxide 
with  stainless  steel  cladding,  only  the  weight  of  the  uranium 
dioxide  is  considered  in  the  unit  MWD/T  (Kef  19:260). 

Dragnev  and  Beets  examine  the  feasibility  of  using 
non-destructive  gamma  spectrometry  measurements  of  irradiated 
fuel  assemblies  to  determine  the  type  of  fuel  (pre-irradiation 
enrichment),  the  history  of  its  irradiation,  and  the 


operational  history  of  the  reactor  during  the  period  when 
the  measured  fuel  was  in  the  core  of  the  reactor  (Ref  20). 
The  irradiated  fuel  is  from  the  German  Vak  reactor;  the  fuel 
elements  are  made  of  pellets  with  1.27  cm  outer  diameter 
and  1.59  cm  length.  A fuel  rod  is  composed  of  two  48-pellet 
sections;  the  rod  diameter  is  1.45  cm.  The  fuel  consists 
of  2.33%  enriched  uranium  in  the  form  of  sintered  uranium 
dioxide;  cladding  is  0.85  mm  Zircaloy-2. 

Since  it  was  planned  to  measure  each  fuel  assembly 
in  only  one  spot  along  its  axis,  measurements  of  the  axial 
distribution  of  burnup  are  made  to  determine  the  best  spot 


for  measurement 


The  ratios  between  the  intensities  of  the 
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605KeV  gamma  peak  (Csl3^)  versus  the  663  KeV  peak  (Csl37) 

f ' ' 

and  the  796  KeV  peak  (Csl3^)  versus  the  663  KeV  peak  (Csl37) 
are  used  as  burnup  monitors.  The  results  can  be  summarized 

/ 

V 

as  follows: 

(1)  The  axial  distributions  of  the  activity 
ratios  follow  qualitatively  the  distribution 
of  the  integrated  neutron  flux  at  different 
points  of  the  reactor  core. 

(2)  The  axial  distribution  of  Csl3^  activities 
decreases  steeply  to  the  ends  of  the  fuel 
assemblies;  this  rapid  decrease  is  caused  by 
the  fact  that  Csl3^  accumulation  is  nearly 
proportional  to  the  second  degree  of  the  integrated 

( neutron  flux. 

(3)  The  activity  ratio  0813*1/03137  as  a burnup 
monitor  gives  the  correct  burnup  value  despite 
the  lack  of  fuel  in  the  center  of  the  rod 
(void  between  the  two  sections  described  above) 

Fig.  2 shows  the  results  of  a single  point  measure- 
ment of  the  gamma  spectrum  of  a Vak  fuel  assembly  after 
1.58  years  cooling  time.  The  peak  labeled  605  is  the 
605  KeV  peak  of  Csl3*l;  the  other  peaks  of  Interest  are 
662  (Csl37)  and  796  (CS13*l,  but  mislabeled  as  Csl37). 

The  quality  of  the  figure  is  poor  because  it  was  duplicated 
from  microfiche.  The  areas  under  the  full-energy  peaks 
are  related  to  the  amount  of  corresponding  gamma  radioactive 
fission  products  in  the  measured  fuel  assembly. 


K C ( 


Single  point  measurements  of  the  activity  ratios 
Csl34(605  KeV)/Csl37  and  Cslj4  (796  KeV)/Csl37 
for  eleven  Vak  fuel  assemblies  with  cooling  time  of  0.42 
years  are  reported.  Fig.  3 plots  these  measurements 
against  the  calculated  burnup  of  each  fuel  assembly;  the 
circled  points  represent  the  reported  data.  As  the  reported 
data  represents  a limited  range  of  burnups,  a linear  least- 
squares  fit  of  the  data  was  calculated  (Ref  21:43).  The 
points  marked  by  x represent  a projection  of  this  fit.  The 
data  associated  with  the  605  KeV  peak  has  a linear  correla- 
tion coefficient  (R)  of  0.737,  while  for  the  796  KeV  data, 

R=0 . 899 . 

Both  the  initial  and  final  amounts  of  U235  for  each 
fuel  assembly  are  reported,  permitting  calculation  of  U235 
depletion.  Fig.  4 plots  U235  depletion  against  the  two 
activity  ratios.  Again,  circled  points  represent  reported 
data,  and  x's  represent  a least-squares  fit.  For  the  605 
KeV  data,  R=0.706;  for  the  796  KeV  data,  R=0.896. 

Using  the  796  KeV  data,  Dragnev  and  Beets  formulated 
a method  for  making  the  Csl34/Csl37  ratio  appear  to  be  a 
far  better  indicator  of  burnup  than  demonstrated  by  Fig.  3. 
Their  treatment  Involves  dividing  each  activity  ratio  by 
the  activity  ratio  of  the  furl  assembly  with  highest  burnup, 
and  dividing  each  burnup  by  * he  highest  burnup;  they  calculate 
a linear  correlation  factor  of  0.993  for  the  relationship 
between  relative  burnup  and  relative  activity  ratio. 

Bresesti  and  Peronl  report  on  experimental  work  corre- 
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605  KeV  (Cs134)/662  KeV  (Csl37) 


Csl3VCsl37  ACTIVITY  RATIO  x 10 


Fig.  3.  Burnup  Versus  Csl3VCsl37  Activity  Ratio  at  0.42  Years 
Cooling  Time  for  Different  Gamma  Lines  of  Csl34 


Csl3VCsl37  ACTIVITY  RATIO 


U235  depletion  versus  Csl3VCsl37  activity  ratio  at 
years  cooling  time  for  different  gamma  lines  of  Csl34 
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lating  the  0813^/03137  activity  ratio  to  the  ratio  Pu/U  in 
irradiated  fuel  assemblies  (Ref  23);  fortunately,  their  re- 
sults are  presented  In  a manner  which  permits  the  Csl3VCsl37 
activity  ratio  to  be  correlated  to  burnup.  They  measure 
the  gamma  spectrum  of  twelve  fuel  assemblies  in  the  spent 
fuel  pond  of  the  Trino  PWR  reactor;  cooling  time  is  not 
reported.  The  burnup  values  are  generated  by  a two-dimen- 
sional burnup  code  (BURSQUID)  using  three  neutron  energy 
groups  (2  thermal,  1 fast);  previous  experiments  validate 
the  accuracy  of  the  code.  While  Dragnev  and  Beets  measure 
each  assembly  at  only  a single  point,  Bresesti  and  Peroni 
measure  each  corner  of  the  assembly  at  nine  levels.  The 
countings  for  the  nine  points  are  summed,  and  then  the 
areas  of  the  photopeaks  for  the  four  corners  are  summed 

to  give  an  integral  value  for  the  entire  fuel  assembly. 

» 

This  integral  spectrum  for  the  fuel  assembly  corresponds 
to  a total  counting  time  of  96  minutes. 

Three  sets  of  activity  ratios  are  reported:  Csl3^(605 
KeV)/Csl37  (663KeV),  Csl3Jl  (796  + 802  )/Csl37  ( 662 ) , Csl3^ 

( 1365 )/CS137 ( 662 ) . These  ratios  are  plotted  versus  burnup 
in  Fig.  5;  the  circles  represent  data  from  fuel  assemblies 
with  an  initial  enrichment  of  3.13/S,  and  the  x's  represent 
fuel  assemblies  with  an  Initial  enrichment  of  3.905?.  The 
lines  represent  a least  squares  fit  to  the  data;  R values 
are:  0.8393  for  the  605  KeV  data,  0.5679  for  the  (796+802) 

KeV  data,  0.69^6  for  the  136*3  KeV  data. 

Bresesti  and  Peroni  also  correlate  Csl37  activity  to 
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Csl3VCsl37  ACTIVITY  RATIO 


Fig.  5.  Burnup  Versus  Csl3VCsl37  Activity  Ratio  for  different 
Gamma  Lines  of  Csl3^ 
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burnup;  R is  0.865  for  this  correlation.  This  relationship 
is  of  interest  because  the  next  experiment  to  be  examined 
deals  with  correlation  of  the  activity  ratio  Csl34/Csl37 
to  Csl37  activity. 

Leender  assumes  a linear  relationship  between  burnup 
and  Csl37  activity;  his  experiment  examines  the  relation- 
ship between  Csl37  activity  and  the  activity  ratio  Csl34/ 
Csl37  (Ref  24).  He  measures  the  gamma  spectrum  of  two 
sets  of  fuel  rods  from  the  Sena  reactor;  one  set  has  half 
the  burnup  of  the  other  (actual  burnups  not  specified). 

Cs134  activity  is  considered  to  be  given  by  the  mean  of  the 
605  and  796  KeV  peaks.  The  relationship  between  Csl37 
activity  and  the  Csl34/Csl37  activity  ratio  for  fuel  rods 
with  the  lesser  burnup  has  an  R of  0.9934.  R for  the  greater 
burnup  rods  is  0.9955. 

Beets,  Bemelmans  and  Pirard  examine  the  Csl34/Csl37 
activity  ratio  correlation  to  burnup  in  two  phases  (Ref  25). 
The  first  phase  involves  the  non-destructive  gamma  spectro- 
metric  measurement  of  fuei  assemblies  in  the  spent  fuel 
pool,  while  the  second  phase  deals  with  destructive  measure- 
ment (fuel  is  decladded  and  dissolved).  This  experiment 
was  intended  to  determine  the  effect  of  mutual  shielding 
and  attenuation  when  entire  fuel  assemblies  are  measured. 

In  the  non-destructive  phase,  the  gamma  spectrum  for  each 
fuel  assembly  is  measured  in  a single  point,  at  the  middle 
of  the  active  length;  the  time  of  measurement  is  15  minutes. 
To  account  for  attenuation  and  detector  efficiency  effects, 
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a correction  factor  F is  determined  for  each  measurement. 
F is  given  by  the  relation 


F=  of 3 (i G-nt  + deoi)/ Cl 


where  the  C's  are  the  areas  under  the  Csl3^  peaks.  A 
linear  correction  is  then  applied  to  the  measured  Csl37 
activity  (662  KeV  peak)  using  the  relation 

C=  F[cCiU  - c^)/(cm  / - c^JaZs  < 4 > 

where  A is  activity,  and  C662  is  the  area  under  the  662  KeV 
peak  of  Csl37-  The  activity  ratio  Csl3VCsl37  refers  to 
the  measured  Csl3^  activity,  which  is  determined  using  the 
605  KeV  peak,  and  the  calculated  Csl37  activity,  which 
is  defined  by  Eq(^).  This  activity  ratio  is  correlated  to 
the  weighted  burnup,  where  weighted  burnup  is  defined  by 


the  expression: 


BVJ=  £ b(vJl/\aI ) e 


0 > 


where 


BW  = weighted  burnup 

?t.  = number  of  months  of  fuel  cycle 
B = burnup 

IaJl  = monthly  integrated  power 
I/s/  = total  integrated  power 
ti  = time  of  irradiation 
7 L.  = time  of  measurement 

fn 
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The  non-destructive  phase  of  this  experiment  results  in  an 
R of  0.9946  for  the  relationship  between  the  Csl34/Csl37 
activity  ratio  and  the  weighted  burnup. 

The  destructive  phase  involves  chemical  decladding 
and  dissolution  of  the  fuel  with  nitric  acid;  the  gamma 
spectrum  of  the  solution  is  then  measured.  The  activity 
ratio  and  weighted  burnup  are  detennined  in  the  same  manner 
as  in  the  non-destructive  phase;  R for  the  destructive 
measurement  phase  is  0.9894. 

Gualandi  presents  interesting  results  on  the  corre- 
lation between  the  Csl34/Csi37  activity  ratio  and  the 
mass  ratio  Pu/U  (Ref  26).  Unfortunately,  his  reported 
data  is  too  incomplete  to  permit  Its  application  to  a burnup 
correlation  as  was  done  with  the  work  of  Dresesti.  However, 
reexamination  of  the  Bresesti  data  (Ref  23-’4)  demonstrates 
a linear  relationship  between  burnup  and  the  Pu/U  ratio 
(R=0.9992) ; this  justifies  qualitative  application  of 
Gualandi's  results  to  the  correlation  between  the  Csl34/Csl37 
activity  ratio  and  burnup.  Gualandi  shows  that  initial  fuel 
enrichment  has  only  a slight  influence  on  the  Csl34/0sl37 
and  Pu/U  correlation  (Rig.  6),  while  power  history  has 
a strong  influence  (Fig.  7).  He  concludes  that  safeguards 
application  (on-site)  of  the  Csl34/CS137  activity  ratio 
requires  a thorough  knowledge  of  power  history.  Since  the 
absence  of  Information  about  power  history  is  the  premise 
upon  which  the  requirement  for  an  off-site  monitor  is  based, 
the  degree  of  power  history  dependence  Is  of  fundamental 
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Pig.  7.  Influence  of  the  Power  Rating  on  the  Correlation  between 
Cs134/Cs137  and  Pu/U  (Ref.  26:619) 


importance.  The  power  history  dependence  of  the  Csl34/Csl37 
activity  ratio  as  an  indicator  of  burriup  is  analyzed  in 
the  following  chapter. 

Further  information  on  the  on-site  application  of 
the  Csl3^/Csl37  activity  ratio  is  available  in  references 
27  and  28;  Ursu  examines  this  ratio  as  an  indicator  of 
cooling  time  (Ref  29). 

The  literature  survey  yielded  only  one  experiment 
dealing  with  a neodymium  ratio  as  a burriup  monitor;  fortun- 
ately it  examined  the  one  ratio  of  Interest  herein,  Ndl46/Ndl45. 
Before  presenting  an  analysis  of  that  experiment,  it  is 
necessary  to  define  a third  measure  of  burriup,  atom  percent 
fission.  Atom  percent  fission  (Ft)  is  defined  as  the  number 
of  fissions  divided  by  the  initial  number  of  heavy  atoms 
times  100(Ref  16:164).  Determination  of  Ft  burriup  requires 
measurement  of  the  number  of  atoms  of  a fission  product 
monitor  and  the  number  of  residual  heavy  atoms;  the  fuel 
specimen  is  dissolved  and  mass  spectrometry  used  for  measure- 
ment of  number  densities.  Burnup  Is  then  calculated  from 
the  relationship: 

Ft  = 100  /CM/y)/(H+M/y)J  (6) 

where 

= number  of  atoms  of  fission  product  monitor 
y = fractional  fission  yield  of  M 
l~/  = number  of  residual  heavy  atoms 
The  most  exact  fission  product  monitor  is  Ndl4  8 (Ref  22:59). 
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Koch  examines  the  correlation  between  Ft  and  the  mass 
ratio  Ndl46/Ndl1t5  (Ref  30);  he  uses  Ndl48  as  the  fission 
product  monitor  in  determining  Ft.  Koch  calculates  a least- 
squares  fit  for*  measurements  made  on  Trino  fuel  samples 
and  reports  the  equation;  his  results  are  shown  in  Fig.  8. 
The  circles  represent  data  for  the  Trino  fuel  samples;  the 
x's  represent  data  for  Sena  fuel  samples.  The  Sena  fuel 
data  is  not  used  in  the  linear  regression. 

In  summary,  the  linear  relationship  between  the  Csl34/ 
Csl37  activity  ratio  and  U235  depletion  is  well  documented 
experimentally,  although  many  measures  of  U235  depletion 
are  used.  Experimental  validation  of  a linear  relationship 
between  the  mass  ratio  Ndl46/Ndl45  and  U235  depletion  is 
sparse  but  convincing. 

Gualandi  (Ref  2 6)  raises  the  issue  of  a power 
history  dependence  when  examining  the  value  of  the  cesium 
ratio  as  a burnup  indicator;  power  history  dependence  is 
best  defined  by  using  fluence  as  a measure  of  burnup. 

For  constant  flux,  fluence  is  the  product  of  flux  and  time; 
a particular  fluence  can  be  obtained  by  different  combina- 
tions of  flux  arid  time.  Gualandi  contends  that  the  Csl34/ 
Csl37  activity  ratio  Is  a function  of  flux  and  time  rather 
than  an  independent  function  of  fluence.  The  degree  of 
this  power  history  dependence  Is  critical  to  the  utility 
of  the  cesium  ratio  as  an  off-site  indicator  of  burnup. 
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IV.  Power  History  Dependence  of  the  Activity  Ratio 
Csl3^/Csl37  as  an  Indicator  of  Burnup 


Determination  of  the  degree  of  power  history  dependence 
of  the  activity  ratio  Csl3VCsl37  as  an  indicator  of  burnup 
requires  a computer  code  which  can  calculate  the  accumula- 
tion arid  decay  of  fission  products  for  varying  power  histories. 
The  accumulation  and  decay  of  fission  products  in  fuel 
elements  can  be  described  by  a set  of  linear  differential 
equations  of  the  form  (Ref  18:19): 


d A lj  jdt  — dL  Nji  Of  & j +-  A (j  t <%  <t>  -f-  A A/j 

- A l- 01  <f>  - XN- 


(7) 


where 


. = number  of  atoms  of  isotope  of  atomic  nujnber 
J i and  mass  J 


= sum  over  all  fissile  isotopes  present  in  fuel 


N* 


it  - 

t 

CTf- 


number  of  atoms  of  fissile  isotope  of  atomic 
number  k and  mass  1 


fission  cross  section  (barns) 

direct  fission  yield  of  fission  product  (i,j) 
neutron  flux  (cm”  sec  ) 
thermal  capture  cross  section  (barns) 

^ = decay  constant  (sec-1) 

For  each  fission  product  isotope  a separate  differential 
equation  exists  and  these  equations  are  coupled  via  and 


*r 

4>  = 
a ; - 
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A/j  ; these  differential  equations  can  be  solved  analytically 
by  Laplace  Transform  for  a time  constant  0. 

Several  assumptions  are  required  to  develop  a simple 
fission  product  inventory  code;  the  first  of  these  assumptions 
is  that  the  neutron  flux  is  time  constant.  Actually,  the 
irradiation  history  In  a reactor  consists  of  a series  of 
constant  power  cycles  with  the  neutron  flux  varying  as  a 
function  of  the  fissile  inventory.  The  second  necessary 
assumption  is  that  U235  is  the  only  fissile  isotope  in 
the  fuel;  this  assumption  disregards  fast  fission  of  U238 
and  also  the  fission  of  accumulating  plutonium. 

Fig.  9 shows  the  production  scheme  for  Csl3iKRef  31)* 
Tables  I and  II  show  the  nuclear  data  required;  Table  III 
defines  the  symbols  used  In  the  two  preceding  tables 
(Ref  2:65,72,79). 

The  final  assumption  is  that  capture  cross  sections 
less  than  one  barn  can  be  disregarded;  the  validity  of  this 
assumption  will  be  examined  when  the  code  is  completed. 

The  largest  capture  cross  section  in  the  mass  132  chain  (Fig. 
9)  Is  that  of  Xel32;  since  this  cross  section  is  less 
than  1 barn,  only  the  mass  133  chain  is  required  in  the 
production  scheme  of  Csl3^«  The  coupled  equations  which 
must  be  solved  are 


d^/33  / ^ - ^tiS  ^ X/33  0 crt  0 / X N /33 


a Ay 


si, 

/S3 


(8) 


30 


Cs  134 


c,-m 


Sb-1 33 


Fig.  9 • Production  Scheme  for  CslB1* 


Nuclear  Data  for  Fissile/Isotopes  (Ref.  2:72) 
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Nuclear  Data  for  Cesium  Isotopes  and  their  Precursors 
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where  0"^  is  the  absorption  cross  section.  These  equations 
are  solved  in  Appendix  B.  These  tedious  hand  calculations 
can  be  avoided  by  using  the  MIT  code  MACSYMA  2 72;  this 
code's  solution  of  Eq  (14)  is  shown  in  Pig.  10. 

The  production  scheme  for  Csl37  is  shown  in  Pig.  11 
(Ref  31)-  Since  the  capture  cross  sections  for  the  members 
of  the  mass  136  chain  are  all  less  than  1 barn,  only  the 
mass  136  chain  is  required  in  the  production  scheme  of 
Csl37.  In  addition  to  Eq  (14),  the  coupled  equations  which 
must  be  solved  are: 

= N0iis  a4  4 + A/,*  CTc  <f>  + AN^ 

- A/,f,  Crc  (f>  - AM,f7  05) 

~ (r,37  <t>  4-  A/,Jt  Crc  <f>  4-  AA/(J7 

- fd,i?  4>  ~ A (16) 

d^Jdi  = Nzk  crf  & <£>  4 hl,ft  <Z  $ f A A/** 

- /d,f7  crcf  - AA/;f)  a?) 

The  solution  to  these  equations  is  contained  in  Appendix  B. 
The  instantaneous  number  of  Csl34  and  Csl37  atoms  are  given 
by  the  expressions 

,Ci  /7  —Cxf-  * -Cxt  si  _ ~Oo  t ~ 

A//j^  *"  e C + C-l o ^ -L  Clj.  C. 

+C*e<«t4  Cne^*t  ua) 
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Pig.  10  MACSYMA  272  Solution  of  a Differential  Equation 


(19) 


= ft 


2t 


+ h 


'IS 


-toft  X 

e / b 


13 


e-^ 

e'^ 


where  the  C's  and  D's  represent  constants  defined  In  Appendix 
B,  and  t represents  time. 

Using  these  expressions  for  the  number  of  Csl3^  and 
Csl37  atoms,  a cesium  inventory  code  (TEST)  was  written 
which  calculates  the  amount  or  activity  of  each  of  these 
isotopes  with  input  parameters  of  Initial  amount  of  U235» 
neutron  flux  and  time.  The  first  application  of  the 
code  was  to  compare  it  to  ORIGEN  (Ref  2);  ORIGEN  Is  the 
model  used  for  PWR's  and  BWR's  by  the  U.  S.  Nuclear  Regulatory 
Commission  (Ref3:28).  Tables  TV  and  V are  output  from  ORIGEN 
being  used  to  model  a PWR  running  at  a constant  average 
specific  power  of  30  MW  per  metric  ton  of  uranium  charged 
to  the  core.  The  same  input  parameters  used  by  ORIGEN 
were  used  in  TEST;  the  listing  and  output  are  contained  in 
Appendix  C.  The  output  of  both  programs  is  for  110  day 
increments.  The  depletion  of  U235  shown  by  TEST  Is  within 
1%  of  the  depletion  of  U235  given  by  ORIGEN  for  the  entire 
fuel  cycle  (1100  days).  Fig.  12  compares  the  amount  of 
Csl3^  calculated  hy  each  code.  TEST'S  result  is  36$ 
higher  than  ORTGEN's  at  J 10  days ; the  difference  decreases 
to  13$  at  1100  days.  Since  all  the  assumptions  made  to 
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TABLE  IV 

Isotopic  Concentrations  (G-Atom/Metricton  of  U)  of  Actinides 
and  Their  Daughters  as  a function  of  Irradiation  Time  (Days) 
TIef  2:106) 
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Csl3^  Buildup  Versus  Time  for  Codes  TEST  and  ORIGIN 
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simplify  TEST  tend  to  decrease  the  production  rate  of 
Csl34,  TEST  was  expected  to  calculate  a lower  value  of 
Csl34  than  ORIGEN.  That  TEST  calculates  a higher  value 
than  ORIGEN  causes  the  accuracy  of  the  mathematical 
methods  used  In  ORIGEN  to  be  suspect.  The  relative  decrease 
in  amount  of  Csl34  (+36%  to  +13%)  with  Increasing  time 
is  as  expected.  ORIGEN  accounts  for  the  fission  of 
accumulating  Pu239>  while  TEST  considers  only  U235  fission. 

Pig.  13  compares  the  amount  of  Csl37  calculated  by 
each  code.  TEST'S  result  is  35%  higher  than  ORIGEN '3 
at  110  days  and  4%  lower  at  1100  days.  Again,  the  relative 
decrease  in  amount  of  Csl37  (+35%  to  -4%)  is  due  to  the 
fact  that  TEST  disregards  the  fission  of  accumulating 
Pu239. 

Fig.  14  compares  the  Csl34/Csl37  mass  ratios  calcu- 
lated by  TEST  and  ORIGEN;  there  is  less  than  0.1%  differ- 
ence at  110  days  and  16%  difference  at  1100  days.  The 
comparison  between  TEST  and  ORIGEN  can  be  summarized  as 
follows.  There  Is  close  agreement  on  U235  depletion  and 
Csl 34/Csl37  mass  ratio.  TEST  is  35%  higher  than  ORIGEN 
for  both  isotopes  at  110  days,  but  the  ratio  (TEST  value/ 
ORIGEN  value)  decreases  with  time  for  both  Isotopes.  It 
Is  concluded  that  TEST  is  an  acceptable  code  for  examining 
the  power  history  dependence  of  the  cesium  ratio  as  an 
indicator  of  burnup. 

The  results  of  the  comparison  between  TEST  and  ORIGEN 
can  also  be  used  to  examine  the  impact  of  disregarding  capture 
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cross  sections  of  less  than  one  barn  In  the  design  of  the 
code  TEST.  Appendix  C shows  that  TEST  calculates  8.08E5  grams 
of  Csl37  after  110  days  and  5.00E6  grains  after  1100  days. 

Eq  (17),  which  is  the  production  rate  equation  for  Csl37, 
has  two  loss  terms.  One  loss  term  is  for  capture,  and  the 
other  is  for  decay.  TEST  disregards  the  capture  term 
because  the  capture  cross  section  of  Csl37  is  O.I69  barns. 
Using  the  above  values  for  number  of  Csl37  atoms  and  Eq  (17), 
loss  rates  due  to  capture  and  delay  were  calculated.  At 
110  days,  the  capture  loss  rate  is  3. 52E-6(grams/sec ) and 
the  decay  rate  is  5-92E-4;  at  1100  days,  the  capture  loss 
rate  is  2.18E-5  and  the  decay  rate  is  3.66E-3-  In  both 
cases,  the  capture  loss  rate  is  two  orders  of  magnitude 
less  than  the  decay  rate. 

All  applications  of  the  code  TEST  hereafter  described 
use  nuclear  data  from  Tables  I and  II  and  an  initial  U235 
charge  of  1.40E2  metric  tons  (same  as  ORIGEN  PWR  model). 

The  first  application  of  TEST  to  the  question  of  power 
history  dependence  Involved  constant  flux  runs  to  a fluence 
of  8.94  E20  (neutrons  per  square  cm).  Nine  flux  levels 
(from  ORIGEN  PWR  model  operating  at  30  MW  per  metric  ton 
of  U)  were  used;  irradiation  time  required  for  the  desired 
fluence  was  300  days  for  the  highest  flux  and  400  days  for 
the  lowest  flux.  Table  VI  lists  the  flux,  time,  and 
activity  ratio  Csl34/Csl37  for  a fluence  of  8.94  E20;  a 
25#  change  In  flux  shows  only  a 1.5%  change  in  activity 
ratio. 
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TABLE  VI 


Activity  Ratio  Csl3VCsl37  (R)  for  Different  Power  Histories 
to  a Fluence  (phi  x t)  of  8.9^  E20  (phi  in  neutrons  per  square 
cm  sec;  t in  sec) 
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TABLE  VII 

Activity  Ratio  Csl3VCsl37  (R)  for  Different  Power  Histories 
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A second  application  of  TEST  examined  more  extreme 
variations  in  flux  level;  a listing  and  the  output  of  this 
program  are  contained  in  Appendix  D.  The  purpose  of  this 
program  was  two-fold;  only  the  odd  lines  of  output  are 
addressed  at  this  point.  Ten  flux  levels,  from  1E13  to 
1E14  (neutrons  per  square  cm  sec),  were  run  for  ten  time 
periods  each;  the  time  periods  ran  from  30  to  300  days,  in 
increments  of  30  days.  The  output  lists  the  flux,  time 
(sec),  and  activity  ratio  Csl  3VCsl37 . Recalling  the 
theoretical  correlation  between  fluence  and  Csl34/Csl37 
activity  ratio  addressed  In  Chapter  IT,  It  is  interesting  to 
note  that  TEST  shows  a linear  correlation  coefficient  of 
0.9999  for  this  relationship  for  a constant  flux.  Returning 
to  power  history  dependence.  Appendix  D shows  that  for  a 
fluence  of  2.59E20,  changing  the  flux  by  a factor  of  2 
changes  the  activity  ratio  by  0.7$.  For  the  same  fluence, 
changing  the  flux  by  a factor  of  S changes  the  activity  ratio 
by  12$,  and  a factor  of  10  changes  the  activity  by  23$. 

In  other  words,  if  the  Csl34/Csl37  activity  ratio  Is  used 

as  a monitor  of  fluence  and  the  uncertainty  in  flux  level  is  a 

factor  of  5,  then  the  uncertainty  in  fluence  will  be  12$. 

Another  way  in  which  power  history  can  vary  is  the  se- 
quence and  duration  of  power  shutdowns  in  fuel  cycle. 

Table  VII  shows  the  TEST  output  for  two  power  cycles  with 
terminal  fluence  of  0.91E20;  the  flux  levels  are  those 
used  by  the  0RIGEN  PWR  model.  Case  1 has  a 30-day  shutdown 
in  the  fourth  month,  while  Case  2 has  a 30-day  shutdown 
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in  the  seventh  month.  The  output  lists  the  activity 
f ratio  at  the  end  of  each  month;  the  difference  in  activity 

ratios  at  the  terminal  fluence  is  0.7%.  A similar  program 
was  run  to  determine  the  effect  of  shutdowns  of  different 
lengths;  two  fuel  cycles  of  constant  flux  and  3G0-day  dura- 
tion were  modeled.  The  first  case  had  a flux  of  7E13 
with  a 30-day  shutdown  in  the  fourth  month;  terminal 
fluence  and  activity  ratio  for  case  1 were  1.633E21  and  1.97, 
respectively.  Case  2 used  a flux  of  8E13  with  a 60-day 
shutdown  in  the  fourth  and  fifth  months;  terminal  fluence 
and  activity  ratio  were  1.659E21  and  1.99,  respectively. 

The  greater  fluence  (1.6/5)  has  an  activity  ratio  1%  higher. 

Tiie  power  history  dependence  of  the  0313^/05137 
activity  ratio  as  an  Indicator  of  fluence  can  be  summarized 
as  follows.  Shutdowns  whose  length  is  short  compared  to 
the  half-life  of  0313^  (2  yrs.)  have  an  insignificant 
impact  on  the  accuracy  of  the  cesium  ratio  as  an  indicator 
of  fluence.  However,  an  order  of  magnitude  uncertainty 
in  the  flux  level  of  the  reactor  results  in  a 23%  uncertainty 
in  the  fluence  Indicated  by  a particular  cesium  ratio. 

As  mentioned  earlier,  the  program  contained  in 
Appendix  D was  designed  to  serve  a dual  purpose.  The  odd 
lines  of  output  are  for  a Csl3^  direct  fission  yield  of 
zero,  and  the  even  lines  of  output  are  for  a CS131*  direct 
fission  yield  of  5.2E-6.  Goodwin  cites  the  5*2E-6  value 
(Ref  3:112);  ORIGEN  uses  zero.  The  difference  in  amount 
of  Cal3^  calculated  by  TEST  using  these  two  values  Is 
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significant  for  low  fluences;  a flux  of  1E13  for  30  days 
results  in  a difference  in  Csl34  calculated.  An  order 
of  magnitude  increase  in  fluence  decreases  the  difference 
l ln  CslB^  calculated  by  an  order  of  magnitude. 


V.  Escape  from  Fuel  to  Environment 

f 

Analysis  of  the  escape  from  fuel  to  environment  of 
the  isotopes  of  interest,  CslB^/lB?  and  Ndl^5/1^6,  is 

V 

addressed  in  three  parts:  escape  from  fuel  matrix,  breach 
of  fuel  cladding,  escape  from  coolant  to  environment. 

The  rate  at  which  fission  products  are  able  to  breach  the 
cladding  is  determined  by  the  burnup  of  the  fuel;  since 
burnup  detection  is  the  ultimate  goal,  breach  of  cladding 
will  be  examined  most  thoroughly. 

Escape  from  Fuel  Matrix 

Pig.  15  shows  the  possible  fission  product  escape 
mechanisms  from  the  fuel  matrix  to  the  fuel  rod  plenum 
^ (Ref  32:5).  The  fuel  matrix  is  a ceramic,  sintered  uranium 

dioxide;  this  ceramic  is  formed  into  pellets  and  stacked 
inside  a Zircaloy  tube.  A helium-filled  plenum  is  provided 
at  the  top  of  the  tube  (rod)  to  collect  fission  product 
gases.  The  ends  of  the  tube  are  sealed  with  welded  plugs 
(Ref  33:0).  Prior  to  sealing,  a helical  compression 
spring  is  placed  to  hear  on  top  of  the  pellet  stack  to 
prevent  shifting  of  f ue L pellets  during  handling  (Ref  3^:^?). 
A temperature  as  high  as  1800  degress  C can  exist  in  the 
center  of  a new  fuel  pellet  operating  at  a linear  heat 
generation  rate  of  13-^  kW/ft(Ref  35:310),  while  the 
primary  coolant  in  a typical  LWR  operates  at  300  degrees  C. 


There  are  four  possible  escape  mechanisms  from  fuel 


to  plenum  (Ref  32:7): 

(1)  diffusion,  trapping,  and  release  from  traps 

(2)  knockout  release  and  surface  effects 

(3)  bubble  migration 

(4)  release  through  cracks  in  the  fuel  matrix 

In  the  temperature  range  800  to  3.800  degrees  C, 

fission  product  movement  occurs  by  a complex  diffusion 
trapping  mechanism;  the  fission  product  atom  diffuses 
through  the  fuel  material  until  trapped.  There  are  two 
categories  of  traps.  Intrinsic  traps  and  point  defects. 
Intrinsic  traps  are  flaws,  such  as  grain  boundaries,  in 
the  fuel  material;  they  are  an  inherent  property  of  the 
I material.  Point  defects  are  formed  by  irradiation  and 

annealed  out  by  heat.  The  trapped  fission  product  atom 
is  eventually  released  as  a result  of  thermal  agitation. 
This  process  continues  until  the  fission  product  atom 
reaches  the  outside  of  the  fuel  material.  The  diffusion 
trapping  mechanism  is  the  dominant  escape  mode  because  it 
occurs  In  the  temperature  range  where  LWR  fuels  operate. 

If  a fission  occurs  within  20  Angstroms  (A)  of 
the  surface  of  the  fuel  material,  the  kinetic  energy  of  a 
fission  fragment  may  drive  it  through  the  surface;  some 
of  the  surface  molecules  are  knocked  out  and  trapped 
fission  products  released.  This  escape  mechanism  is  called 
knockout  release;  below  700  degrees  C,  knockout  release 
is  the  dominant  escape  mode. 
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Bubble  migration  occurs  at  high  fuel  temperatures 
(above  1800  degrees)  when  the  fuel  matrix  begins  to  lose 
strength;  bubbles  form  when  diffusing  atoms  Joint  together, 
v Bubble  migration  occurs  along  a thermal  gradient  by  a 

mechanism  called  surface  diffusion.  The  formation  and 
movement  of  these  bubbles  causes  swelling  in  the  fuel 
matrix.  Bubble  migration  is  not  a significant  escape 
mode  because  LWR  fuels  do  not  reach  a high  enough  temperature. 

Rapid  temperature  changes  in  the  fuel  matrix  can 
cause  the  low  tensile  strength  ceramic  to  crack;  diffusing 
fission  products  can  then  be  released  to  the  plenum  by 
these  cracks.  Since  rapid  reactor  start-up  or  shutdown  is 
not  a normal  practice,  release  through  cracks  is  not  a 
significant  escape  mechanism  during  routine  reactor  operation. 

Since  diffusion  (including  trapping  and  release) 
is  the  dominant  escape  mode  from  fuel  matrix  to  plenum, 
it  is  appropriate  at  this  point  to  examine  the  diffusivity 
of  cesium,  neodymium  and  their  precursors.  Goodwin's 
assumption  that  the  use  of  a same-element  ratio  ensures 
uniform  fractional  release  from  the  fuel  matrix  (Ref  3:^9) 
is  not  Justified;  the  fractional  release  of  any  member 
of  a decay  chain  is  not  simply  and  exclusively  determined 
by  its  own  diffusivity  (Ref  36:89).  The  fact  that  each 
precursor  is  able  to  diffuse  before  decaying  will  increase 
the  fractional  release  of  a particular  isotope  above  the 
level  which  estimates  ignoring  the  mobility  of  precursors 


1 


% 


would  suggest.  Decay  products  having  precursors  with 
relatively  high  dif fusivities  and  half-lives  which  are  a 
significant  fraction  of  irradiation  time  will  experience 
the  greatest  increase  in  fractional  release. 

The  precursors  of  both  neodymium  Isotopes  of  interest 
are  cerium  (Ce)  and  praseodymium  (Pr);  only  Prl^5  has  a 
half-life  of  significant  length,  6 hours  (Ref  11:98). 

There  is  little  to  be  found  in  the  literature  about  the 
diffusivity  of  non-volatiles;  the  boiling  point  of  neodymium 
and  both  its  precursors  Is  over  3 Q degrees  C.  However, 
Brown  and  Faircloth  state  that  cerium  is  relatively  immobile 
(Ref  37:31),  and  they  conclude  that  cerium  escapes  the 
fuel  matrix  only  as  a recoiling  fission  fragment.  This 
conclusion  implies  that  only  cerium  atoms  produced  with 
10A  of  the  surface  of  the  fuel  matrix  have  a chance  of 
escaping  (Ref  3?:9)  • It  is  logical  to  assume,  based  on 
non-volatility,  that  Pr  and  Nd  are  as  immobile  as  Ce. 

Since  neither  Pr  nor  Nd  Is  a direct  fission  product 
(Ref  2:80),  the  only  Nd  to  be  found  outside  the  fuel 
matrix  is  that  formed  by  decaying  Ce  which  lias  escaped 
by  recoil  as  a fission  fragment.  For  one  of  the  Vak  fuel 
pellets  described  iri  Chapter  III,  only  3-15E-7  of  the 
volume  lies  within  10A  of  the  surface.  Even  if  one  assumes 
that  all  cerium  atoms  produced  in  tills  volume  recoil 
outward,  only  3.15E-7  of  the  cerium  atoms  produced  within 
the  fuel  can  escape  the  fuel  matrix.  From  the  standpoint 
of  escape  rate,  this  effectively  reduces  the  neodymium 
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fission  yield  by  a factor  of  3-15E-7. 

The  volatility  of  cesium  gives  it  great  mobility 
within  the  fuel  matrix;  it  migrates  radially  and  axially 
in  the  direction  of  decreasing  temperature  (Ref  38).  Both 
cesium  isotopes  of  Interest  have  iodine  and  xetion  as 
precursors,  but  the  only  half-lives  of  significant  length 
in  either  chain  are  those  of  1133(21  hr.)  and  Xel33 ( 5 • 27d ) . 
Since  both  iodine,  a halogen,  and  xenon,  a noble  gas,  have 
greater  dif fusivities  than  cesium  (Ref  39),  the  fractional 
release  of  Csl3^  should  be  greater  than  that  of  Csl37; 

Brown  and  Faircloth  investigate  and  verify  this  experi- 
mentally (Ref  37).  Friskney  and  Speight  perform  calcula- 
tions to  quantify  this  difference  in  fractional  release 
(Ref  36:90);  their  results  are  shown  in  Fig.  16.  Curves 
(7)  and  (8)  apply  to  Csl37  because  of  the  short  half-lives 
of  its  precursors.  The  othei*  curves  apply  to  Csl3^.  D1 
represents  the  diffusion  coefficient  of  1133,  and  D2  refers 
to  Xel33-  D3  refers  to  Csl3^,  and  the  unsubscripted  D 
refers  to  Cs.137.  Since  a diffusion  coefficient  is  a function 
of  grain  radius  and  temperature,  two  arbitrary  D's  are 
investigated,  IE-13  and  lE-l^  (square  cm/sec).  For  the 
worst  case,  Di  arid  1)2  both  an  order  of  magnitude  greater 
than  D3,  there  is  a difference  of  27%  in  fractional  release 
of  the  two  isotopes  at  100  days;  this  case  is  represented 
by  curves  (3)  and  (7). 

The  calculations  of  Friskney  and  Speight  involve  a 
simple  diffusion  model;  it  neglects  both  the  trapping-release 
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Fractional  release  vi  time  for  the  final  stable  isotope  (3) 
in  the  chain: 


(for  the  case  of  ml,  1 J,Xc,  mCs,  k,  * 9 255  X 10  *. 
k2-  1.522  x 10'6/scc). 


Fig.  16.  Fractional  Release  Versus  Time  for  the  Final  Isotope 
(3)  in  the  Chain 
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mechanism  involved  in  diffusion  and  the  fact  that  cesium's 
movement  is  actually  a migration  to  lower  temperature. 
However,  their  relative  values  of  fractional  release  should 
be  valid  if  the  correct  absolute  diffusion  coefficients 
are  used.  The  unavailability  of  these  diffusion  coefficients 
prevents  resolution  of  the  issue  of  different  fractional 
release  of  Csl34  and  Csl37. 

In  summary,  only  a very  small  fraction  of  the  neodymium 
(cerium)  atoms  can  escape  from  the  fuel  matrix,  while 
the  cesium  isotopes  have  a high,  but  apparently  different, 
escape  rate. 

Preach  of  Fuel  Cladding 

The  fuel  rod  cladding  acts  as  a barrier  for  fission 
product  transport  and  as  a container  for  the  fission  gas 
pressure  buildup;  it  is  designed  1,0  prevent  fission  products 
from  escaping  to  the  primary  coolant.  Fig.  15  shows  the 
possible  escape  mechanisms  through  the  fuel  cladding. 

Eckart  calculates  diffusion  rate  of  1.3E-4  cm/year  for 
a noble  gas  in  the  fuel  clad  (Ref  32:29);  the  thickness 
of  cladding  on  the  Vak  fuel  described  earlier  is  8.5E-2  cm. 

It  wouid  take  even  the  most  mobile  fission  product  600 
years  to  diffuse  through  the  fuel  clad;  examination  of  the 
fission  product  content  of  a reactor's  primary  coolant 
demonstrates  that  fission  product  release  must  take  place 
through  defects  in  the  fuel  rod  cladding.  Stainless  steel 
was  used  as  fuel  cladding  In  early  I.WR's,  but  only  zirconium 
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alloys  have  been  used  since  1968  (Ref  40). 
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The  causes  of  fuel  failures  (clad  defects)  for  slightly 
enriched  uranium  dioxide  clad  with  Zircaloy  are  reviewed 
by  Robertson  (Ref  40);  he  groups  the  causes  of  fuel  failure 
into  three  categories:  manufacturing  defects,  external 
causes,  fuel-cladding  interactions. 

Manufacturing  defects  include  faulty  components, 
faulty  content  and  faulty  assembly.  Faulty  components 
shows  a very  low  Incidence  because  of  the  simplicity 
of  the  fuel  rod's  design;  of  the  700,000  rods  irradiated 
in  1971,  none  showed  a fabrication  defect  of  any  kind 
(Ref  41:185).  The  most  significant  type  of  faulty  content 
defect  is  contamination  of  the  fuel  by  a hydrogenous 
compound  such  as  water;  hydrogen  interacts  with  the  Zircaloy 
cladding  by  a process  called  hydrlding.  The  formation  of 
zirconium  hydride  makes  the  clad  brittle  and  causes  a 
clad  volume  increase;  these  two  changes  in  the  clad  can 
lead  to  cracking.  The  most  frequent  incidence  of  faulty 
assembly  involved  end-plug  welds;  but  these  are  an  insignifi- 
cant contributor  to  fuel  failure  as  they  are  usually  detected 
by  the  quality  control  measures  of  the  manufacturer. 

External  causes  of  clad  failure  include  the  following: 
fuel  handling,  foreign  bodies,  fuel  overpower,  vibration 
and  fretting,  cladding  corrosion.  In  many  power  reactors 
fuel  management  requires  that  irradiated  fuel  assemblies 
be  repositioned  at  the  end  of  a cycle;  these  assemblies 
are  then  Irradiated  again  in  a subsequent  cycle.  This 
repositioning  subjects  the  cladding  to  mechanical  forces 
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when  it  is  most  susceptible  to  damage,  i.e.  when  embrittled 
by  irradiation  damage  and  some  degree  of  hydriding.  The 
current  design  of  fuel  handling  equipment  has  practically 
eliminated  fuel  handling  as  a cause  of  failure. 

If  a foreign  body  Is  trapped  between  fuel  rods,  it 
may  be  excited  Into  vibration  by  the  flow  of  coolant; 
the  relatively  thin  cladding  is  then  susceptible  to  fretting 
(chafing)  damage.  Start-up  cleanliness  of  the  primary 
coolant  circuit  precludes  this  failure  mechanism. 

Since  in  operation  the  fuel  expands  more  than  the 
cladding,  some  voidage  is  usually  provided  within  the 
fuel  rod  to  allow  for  this  differential  expansion.  The 
magnitude  and  location  of  the  voidage  is  at  least  roughly 
matched  to  the  expected  temperature  distribution  of  the 
fuel,  which  in  turn  depends  on  the  power  developed  in  the 
fuel.  If  the  power  developed  is  considerably  greater 
than  that  for  which  the  fuel  rod  was  designed,  the  cladding 
can  be  overstressed  to  the  point  of  failure. 

In  some  power  reactors,  vibration  of  the  fuel  assemblies 
has  resulted  in  severe  fretting,  and  even  penetration  of 
the  cladding.  In  some  instances  a complete  string  of  fuel 
assemblies  has  vibrated  against  a tube  containing  it,  while 
in  others  a spacer  grid  or  similar  component  has  vibrated 
against  the  fuel  rods. 

Cladding  corrosion  is  assigned  to  external  causes 
since  the  corrosion  resistance  of  Zircaloy  is  adequate 
to  ensure  good  fuel  rod  performance  as  long  as  the  coolant 
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chemistry  is  well  controlled.  Where  corrosion-induced 
fuel  failures  have  been  reported  the  basic  cause  has  been 
a heavy  build-up  of  deposits  on  the  cladding  surface,  re- 
sulting in  the  Zircaloy  corrosion  occuring  at  a much  higher 
temperature  than  expected  on  the  clad  outer  surface. 

Fuel-cladding  interactions  include  the  following 
mechanisms:  rod  elongation,  fuel  densi ficat ion , cladding 
fatigue,  power  ramp  failures.  Rod  elongation  and  fuel 
densification  shorten  the  useful  life  of  a fuel  rod,  but 
do  not  actually  constitute  defect  mechanisms;  neither  of 
these  mechanisms  has  directly  caused  a reported  clad 
rupture.  Cladding  fatigue  also  has  yet  to  be  blamed  for 
any  fuel  failures  in  power  reactors,  but  clad  failure  by 
fatigue  v/ as  demonstrated  by  some  early  steel-clad  experi- 
mental fuel  rods.  Power  ramp  failures,  however,  have  been 
the  direct  cause  of  reported  clad  ruptures. 

When  reactor  power  is  raised  Initially  the  cladding 
is  undamaged  by  irradiation  and  fuel-cladding  clearances 
are  available  to  accommodate  differential  thermal  expansion. 
If  the  power  is  subsequently  increased  from  a low  level  to 
a high  level,  the  cladding  is  somewhat  embrittled  by 
irradiation  damage  while  the  clearances  have  been  at 
least  partly  taken  up  by  fuel  swelling  and  by  the  cladding 
creeping  inward;  most  cladding  will  slowly  deform  under 
the  coolant  pressure  by  irradiation-accelerated  creep.  Thus 
there  exists  the  potential  for  power  ramps  to  crack  the 
cladding;  defect  probability  in  ramping  failures  increases 
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with  power,  power  increase,  burnup,  and  dwell  time  at  high 
power  (Ret  42). 

A fuel-cladding  Interaction  not  discussed  by  Robertson 
is  fission  product  Induced  corrosion  of  cladding;  Wood 
identifies  iodine  corrosion  as  a contributor  to  power  ramp 
failures  (Ref  42:157).  Cesium  has  also  been  associated 
with  cladding  corrosion  (Ref  38:144). 

Of  greater  importance  than  clad  failure  cause  is  clad 
failure  rate.  Multer  summarizes  European  fuel  performance 
(Ref  43).  For  PWR's  with  Westlnghouse  stainless  steel 
Cad  fuel  rods  and  a burnup  of  28,000  MWD/T,  the  average 
failure  rate  has  been  0.02?.  Zircaloy  cladding  has  re- 
duced the  failure  rate  to  0.01?  in  some  cases,  but  for 
PWR's  the  average  failure  rate  remains  close  to  0.02?.  BWR 
fuel  has  shown  an  average  failure  rate  of  0.01?. 

Williamson  and  Proebstle  examine  American  BWR  experience 
and  recent  BWR  fuel  improvements  (Ref  44);  they  conclude 
that  the  newest  General  Electric  BWR  fuel  is  demonstrating 
a failure  rate  of  0.006?,  Kramer  examines  the  American 
experience  with  PWR  fuel  (Ref  45);  he  shows  an  average 
failure  rate  of  0.02?  and  forecasts  no  foreseeable  improve- 
ment. The  lower  failure  rate  in  BWR  fuel  can  be  attributed 
to  the  thicker  cladding  (.025  inch  for  PWR,  .033  inch  for 
BWR)  and  to  the  lower  coolant  pressure  (2000  PSI  for  PWR, 

1000  PSI  for  BWR). 

In  summary,  broach  of  fuel  cladding  is  a result  of 
clad  failure.  Current  clad  failures  are  usually  attributed 


62 


V 


tg  hydriding  and  power  ramps,  and  fission  product  corrosion 
is  recognized  as  a contributor  to  power  ramp  failures.  The 
current  clad  failure  rate  is  0.0  2%  for  PWR ' s and  0.01/5 
L for  BWR's;  a slight  improvement  in  BWR  fuel  is  forecast. 

Data  on  Combustion  Engineering  Inc.  (C-E)  and  Kraftwerf 
Union  (KWU)  fuels  shows  similar  failure  rates  (Ref  46). 

Escape  from  Primary  Coolant  to  Environment 

Fig.  17  shows  fission  product  escape  mechanisms  from 
primary  coolant  to  environment.  Environmental  Protection 
Agency  (EPA)  surveillance  studies  describe  the  radwaste 
systems  designed  to  prevent  this  escape  and  analyze  their 
effectiveness  (Ref  47,  48).  The  relative  abundances  of 
radionuclides  in  waste  effluents  are  functions  of  (Ref  49:3): 
j (1)  their  abundance  in  the  primary  coolant 

(2)  their  respective  half-lives 

(3)  design  of  the  radwaste  treatment  system 

(4)  waste  treatment  practices 

Radwaste  systems  in  use  in  American  LWR's  are  a compro- 

. 

mise  between  cost  and  effectiveness.  The  release  rate  of 
cesium  by  a current  RWR  could  be  reduced  by  a factor  of 
1000  with  a "maximum"  radwaste  system  costing  twice  as 
much  as  a current  radwaste  system  (Ref  50:49,52);  a rad- 
waste system  described  by  the  EPA  as  a "minimum  treat- 
ment" system  would  have  half  the  cost  of  a current  system 
and  a cesium  release  rate  10Q  times  greater.  For  a PWR, 
a "minimum  treatment"  radwaste  system  would  have  one-fourth 
the  cost  of  a current  system  and  a cesium  release 
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rate  100  times  greater  (Ref  50:59,  6l). 

The  maximum  coolant  canal  concentrations  of  Csl3*t 
and  Csl37  measured  by  the  EPA  at  Dresden  I,  a BWR  with 
0.01%  fuel  failure  rate,  are  1.6E-12  Curles/L  and  4.3E-12 
Curles/L,  respectively  (Ref  10).  The  maximum  stack  release 
rate  of  Csl37  measured  at  Dresden  is  5.5E-11  Curies/sec; 
no  measurement  is  reported  for  0313^.  These  measurements 
were  made  Just  before  the  end  of  a fuel  cycle;  some  of  the 
fuel  rods  were  approaching  maximum  rated  burnup. 

The  maximum  coolant  canal  concentrations  of  Csl3^ 
and  Csl37  measured  by  the  EPA  at  Haddam  Neck,  a PWR  with 
0.02%  fuel  failure  rate,  are  4.0E-13  Curies/L  and  3.0E-13 
Curies/L,  respectively  (Ref  12).  The  maximum  stack  release 
rates  of  Csl3^  and  Csl37  measured  at  Haddam  Neck  are 
^.8E-11  Curles/sec  and  3-3E-11  Curles/sec,  respectively. 
Again,  these  measurements  were  made  well  into  the  fuel 
cycle  when  some  fuel  had  attained  high  burnup. 

Of  the  cesium  concent-rat  ions  reported  above,  the  BWR 
has  the  higher  reported  values  for  both  coolant  canal  and 
stack  gas;  the  BWR  values  wiLl  now  be  used  to  determine 
the  limiting  distance  at  which  the  Csl3VCsl37  activity 
ratio  could  be  measured. 
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VI . Detection 

Cesium  escaping  through  the  reactor  stack  is  a very 
fine  particulate;  particulate  samples  are  normally  collected 
on  a filter  medium  with  an  air  pump  and  a flow-measuring 
device.  Gamma-spectrometric  analysis  for  Csl3^  and  Csl37 
requires  collection  of  at  least  1200  cubic  meters  of  air; 
the  minimum  detectable  level  for  both  isotopes  is  IE-14 
Curies/cubic  meter  (Ref  51:10).  Since  the  clandestine 
collection  of  such  a large  volume  of  air  would  require  a 
reasonable  distance  between  reactor  and  collection  point, 
a collection  point  at  1500  m (1  mile)  will  be  assumed. 

Goodwin  (Ref  3:70)  gives  a good  review  of  various 
models  of  atmospheric  dispersion;  EPA  uses  the  following 
generalized  Guassian  diffusion  equation  (Ref  10:57): 


X*  q/i 


TTCTy  iu  e 


(20) 


where 


^ = ground-level  centerline  concentration 
(Curies/cubic  meter) 


Q = release  rate  (Curies/sec) 


0^  - lateral  and  vertical  dispersion  coefficients  (m ) 


= mean  wind  speed  at  height  of  release  (m/sec) 
/y  = effective  release  height  (m) 

Normalized  concentrations  (Xz/c?)  as  a function  of  distance 
from  the  stack  developed  by  Brookhaven  are  used  in  the 


I 

1 


66 


solution  of  this  equation. 


The  maximum  measured  release  rate  of  Csl37  at  Dresden 
is  5-5E-H  Curies/sec.  Csl34  Is  not  reported,  but  can 
be  assumed  to  escape  with  a release  rate  of  the  same  order 
of  magnitude  as  Csl37.  The  EPA  reports  highest  ground-level 
centerline  effluent  concentrations  for  neutral  atmospheric 
stability  and  6 m/sec  wlndspeed  at  release  height  (Ref  10:48); 
the  stack  height  at  Dresden  is  91  m.  Using  the  dispersion 
coefficients  for  these  conditions,  solution  of  Eq(20)  gives 
a Csl37  concentration  of  1.59E-16  Curies/cubic  meter  at  1 
mile  from  the  stack;  this  concentration  is  two  orders  of 
magnitude  less  than  the  EPA's  minimum  detectable  level. 
Clearly,  there  are  other  ways  to  approach  the  detectability 
of  a given  release  rate;  given  the  atmospheric  conditions, 
the  point  of  maximum  ground-level  concentration  can  be 
calculated  and  assumed  to  be  the  collection  point.  However, 
the  preceding  calculation  is  so  biased  in  favor  of  a high 
Csl37  concentration  (particulate  settling  disregarded) 
that  further  treatment  is  unnecessary. 

The  EPA  recommends  a sample  size  of  3*5L  for  gamma- 
spectrometric  analysis  of  cesium  in  water  (Ref  51:10);  the 
minimum  detectable  concentration  is  IE-11  Curies/L.  Clandes- 
tine collection  of  a sample  of  this  size  precludes  use 
of  the  discharge  canal  itself;  dilution  by  the  body  of 
water  receiving  the  discharge  canal  must  be  assumed. 


The  maximum  measured  coolant  canal  concentrations 
at  Dresden  are  1.6E-12  Curies/L  for  Csl34  and  4.3E-12 


Curies/L  for  Csl37.  The  coolant  canal  discharges  Into  a 
river  whose  average  flow  rate  is  11,000  cubic  feet/sec 
(cfs);  the  dilution  rate  in  such  a situation  is  a function 
of  an  unmanageable  number  of  variables.  EPA , however, 
uses  a very  simplistic  approach  to  arrive  at  a dilution 
factor  of  50  (Ref  10:60).  Another  EPA  report  calculates 
a dilution  factor  of  25  for  a river  with  average  flow 
rate  of  18,000  cfs  (Ref  12:67).  The  lower  dilution  factor 
will  be  used  in  the  interest  of  maximizing  the  cesium 
concentration  at  the  collection  point.  Dilution  of  the 
canal  concentrations  by  a factor  of  25  results  in  sample 
concentrations  of  6.4E-14  Curies/L  for  Csl34  and  1.72E-13 
Curies/L  for  Csl37;  these  concentrations  are  two  orders 
of  magnitude  less  than  EPA’s  minimum  detectable  level. 
Again,  tills  calculation  has  been  biased  toward  a high 
cesium  concentration. 

The  fact  that  the  EPA’s  minimum  detectable  level 
(MDL)  for  cesium  In  water  is  an  order  of  magnitude  greater 
than  the  reported  coolant  canal  level  is  explained  by  the 
fact  that  the  MDL  required  for  measuring  a cesium  ratio 
is  an  order  of  magnitude  greater  than  the  MDL  required  for 
simple  measurement  of  cesium.  The  canal  levels  are  Just 
above  the  MDL  for  simple  measurement  of  cesium. 

The  accuracy  of  this  type  of  measurement  is  addressed 
by  Gans  (Ref  52);  he  submitted  identical  samples  containing 
on  the  order  of  IE-9  Curles/L  of  Csl34  and  Csl37  to  28 
laboratories  in  Germany  and  Holland  that  routinely  perform 
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analytical  work  for  the  nuclear  industry.  He  reports  a 
10$  variation  in  the  measurements  reported;  furthermore, 
he  predicts  errors  on  the  order  of  100$  when  measuring 
concentrations  close  to  the  MDL's  used  herein. 
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VII.  Results,  Conclusions  and  Recommendations 

Using  measured  effluent  concentrations  of  CslB^ 
and  Csl37  from  an  American  BWR  with  a characteristic  fuel 
failure  rate  of  0.01%,  calculated  concentrations  at  assumed 
sample  collection  points  are  two  orders  of  magnitude  less 
than  EPA  prescribed  minimum  detectable  levels.  These 
calculations  are  biased  toward  the  high  side,  and  the 
problem  of  different  release  rates  of  the  cesium  isotopes 
from  the  fuel  matrix  is  disregarded.  If  minimum  radwaste 
systems  are  assumed,  the  sample  concentrations  might  approach 
minimum  detectable  level,  but  measurements  at  this  low  level 
of  concentration  are  of  questionable  accuracy.  It  is  con- 
cluded that  the  activity  ratio  Csl34/Csl37  is  not  an  accept- 
able off-site  indicator  of  burnup. 

The  neodymium  isotopes,  Ndl^b  and  Ndl^6,  have  fission 
yields  of  the  same  order  of  magnitude  as  Csl37,  but  only 
3-15E-7  of  neodymium's  precursor  is  able  to  escape  the 
fuel  matrix.  Cieariy,  the  collection  point  concentration 
of  neodymium  will  be  far  less  than  that  of  cesium;  conversion 
of  the  expected  Csl37  concentration  (using  minimum  radwaste 
system)  to  a number  density  and  reduction  of  this  number 
density  by  a factor  of  3.1c-E-7  gives  an  expected  neodymium 
number  density  on  the  order  of  1 atom/L  for  a liquid  sample 
and  10  atoms/L  for  an  air  sample.  Measurement  of  such 
levels  is  not  possible;  the  ratio  Ndlii6/Ndl45  is  not  accept- 
able as  an  off-site  indicator  of  burnup. 
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Radioactive  noble  gases  have  been  successfully 
measured  by  the  EPA  at  1 mile  from  the  stack;  it  is 
recommended  that  radioactive  isotopes  of  Xenon  and  Krypton 
be  examined  as  indicators  of  burnup. 
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Appendix  A 


Versus  Background  Concentrations 


Molecular  Fraction  of  Kr 
Number  Density  of  Air: 
Number  Density  of  Kr: 
Kr84  Natural  Abundance: 
Kr83  Natural  Abundance: 
Kr84  Number  Density: 

Kr83  Number  Density: 


in  Air:  1.0  E-6 

2.548  E19  cm-3 
2.548  E13  cm-3 
0.5690 
0.1155 

1.45  E13  cm-3 
2 94  E12  cm-3 


Molecular  Fraction  of  Xe 
Number  Density  of  Xe: 
Xel31  Natural  Abundance: 
Xel32  Natural  Abundance: 
Xel34  Natural  Abundance: 
Xel36  Natural  Abundance: 
Xel31  Number  Density 


in  Air:  8.0  E-8 

2.04  E12  cm-3 
0.2118 
0.2689 
0.1044 
0.0867 

4.32  Ell  cm-3 


Xel32  Number  Density: 
Xel34  Number  Density: 
Xel36  Number  Density: 


5.49  Ell  cm-3 


2.13  Ell  cm-3 
1.81  Ell  cm-3 

Maximum  Observed  Stack  Gas  Concentration  of  Kr85  at 
Oyster  Creek:  3.2  E-7  yci/cc 

Maximum  Observed  Stack  Gas  Concentration  of  Xel33  at 
Oyster  Creek:  1.9  E-4  pCi/cc 

Jem/iTy  = XhJ 

^Xc./33  ~~  l‘&2.  sec. 

A / = hcn\JiTr A 

^faffs’  = GjvC^ 

A//e/«  " 4.CSEC> 

For  every  100  atoms  of  1^35  thermally  fissioned,  following 
number  of  atoms  created: 
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Kr83  — 

0.48 

Kr84  — 

1.10 

Kr85  - 

1.50 

Xel31  - 

2.64 

Xel32  - 

4.35 

Xel33  - 

0.360 

Xel34  - 

7.60 

Xel36  - 

6.30 

Kr83  and  Kr84  can  be  expected  to  be  present  in  stack  gas  in 
concentration  on  the  order  of  1E6  per  cubic  centimeter. 

Xel31,  132,  133,  134,  136  can  be  expected  to  be  present  in 
stack  gas  in  concentration  on  the  order  of  1E7  per  cubic 
centimeter . 
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Appendix  B 


Cesium  Production  Equations 


<74  * — Bk  i 

c^L  p/J 

sL  pvj*  - A/O)  = 

S/C  - qAL$N?  = A/CO 

=■  a/c°)/ s--C,a 


= //r-4 

A/  - A/O) 


Huf*  A/co  e F 

235' 


8o 


, the  following  solutions  were  determined.  The 


By  similar  calculations 
constants  are  defined  in  Appendix  C. 


/3  3 

A/,« 


-c.it 


~ c/3  G 7 Wf  ^ C/5- 


A/ 


V 


/53 

/33 


c?23  «"4lf -4y  4f  e_<5*S  4*.  e"4** 

^e'4<  -C*  e**+  4,7  e‘**  + <3» 
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/53 
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y 4r 

e <3*> 

•A  e-^  -/■  csl  eT6** 


A C - 4,  e^7  - 4,  * <5>  ^ 


C7,  „ ^<-73 


/sy 


T \-aHee”r+ 


<Sr 


A/,„  * i>7e 
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/37 


7 W ■/-  i)g 
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Appendix  C 
Program  Tes 


PROGRAM  TEST  (INPUT, OUTPUT) 
REAL  11331,11371 
DIMENSIO)  INTEGRA  6) 

DATA  INTEGR/3'(-0),0,2*<-G>/ 
CALL  3Y5TEMC( 115,INTEGK) 
REA0*,U235I,SIGF,SIGA 
.READ*, SOI  331, YSB133,DSB133 
READ*, TE133I , YTE133»OTE133 
READ*, II 3 31, Y £13 3, Cl 133 
READ*,XE133I, XFl 33C , DXE 1 33 
READ*»Cil33I ,CS133C 


REAP*, CS 13 41, YCS134,DCS134,CS134C 
READ*,  H37I,YI137,DI137_  _ 

READ* , Xcl 371 , DXE 1 37 
READ*  , C31  371 , YG5 1 37  , DCS  1 37  , CS  1 37C 


00  101  1=1,10 

READ* , FHI 

ClA=U235r *SIGF*PHI 
C-1=YS3133*C1A 


C2=SIGA*PHI 

C3=DS3133 


C4=S:11 331 

C5=C1/ (03-C2) 
C6=(C1  + (C2*C4_L 
C7=YT  E 1 33*  CIA 


( C 3*C4 ) ) / (C 3-C  2) 


C6=C3“C5 

C9=C3*G6 

C10=0TE133 

Cil=C7+CJ 


Cl2=TE133I  _ 

C13=C 1 1/ ( C 2-C 10 ) *C 1 2 -C9/I C 3 -C 1 0 ) 
C14=C11/ (C10-C2) 

C15=C j/ (C10-C3) 

C16=Y 1 1 3 3 *C1A 
C17=C10*C13 

C18=C10*C14  

ClS=Cio*C15 

C20  = QI  1 33  _ 

C21=C16fC13 

C22=I 1331  

C23=C21/ (C20-C2) 

C 24=01 9/ ( C 20-  03) 

C25=G17/ (C20-C10) 

C26=C22*C17/(010-C20)-C19/(C3-C20) 

C27=XE133C*PH1 

C2fc=DXE 133  _ 

C2G=C27  *023 

C30=C20*C23  _ 

C31  = C20*C24 

C32=C20*C25 

C33=C20*C26 

C34=XE 1331  

C35=C  3 0 / ( C29- C2 ) 

C36=C31/ (C29-C3)  

C37=C32/ ( C29-C10 ) 

C3C =C 33/ ( C29-C20 ) 

82 


P <0G1-  i T £ r>T 


74/74  CPT=1 


FTi'l  4.6*446 


( 


( 


-( 


039=030/ ( C 2-C  2?) -031/ (C3-C29) +C 32/(C10-C 2? >+C 33/(020 -029) +C 34 

C4Q=C3133C*PHI 

041=029*335 

042=026*  £36 

C43=C26*C  37 

044=026*033 

C45=C2c* J39 

C46=CS 1331 

047=341/(040-02) 

C4c=C42/ (040-03) 

C45=C43/(E 40- 310  ) 

050=344/  (C40-320) 

051=045/ (340-329) 

C 52=041/ (C2-C 40) -C42 / ( C 3-C 40 ) + C 43 / ( 0 10 -C 40 > +044/ (C20-C40 ) ♦ 

1 C45/ ( C?  J-C  40  ) +C46 
C53=Y3S134*C1  1 
C54=CS1 330*PHI 
055=031 34C*PH  I 

C56=DCS134  __  

057=054*047 

056=054*043 

059=054*04 J 

060=054*050 

061=054*051 

062=054*052 

063=031341 

064=055+056 

065=053  + 057 

066=065/ (C64- 02) 

C67=-C2+T 

C6f=05c/(Cb4-03) 

C69=-C3*T 

070=059/(064-010) 

C71=-010*T 

C72  = 060/(C  64-020 ) 

C73=-020*l 
074=061/ (064-029) 

C75=-02  3*  T 
076=062/(064-040 ) 

C77=-C4Q*T 

078=065/(0  2-0  64)  -C5b/ ( C3-C64 ) +059/(010-0  64) +060/ (C 20-064)  + 
1061/(029-064) +062/ (040-064)+ 06 3 
C79=-C64*T 

AMT134=C56*EXP(C67)-C63*EXP(C69)+C70*EXP(C71)+C72*EXP{C73)+ 

1074*EXP(075)+076+EXP(C77)+C7b*EXP(C79) 

01=01  A 
D2=C2 
03=YI137 
04=01*03 
05=01137 
06=11371 
07=04/ ( J5-02) 

08=06+04/(02-05) 

09  = 0X  E 1 3 7 
010=05*97 
D11=D5* JJ 
012=XE137I 
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74/74  0PT=1 


i 

pwg^;m 


TEST 


FT N 4.6+446 


0 


013=010/(09-02)  

014=D11/(D9”D5) 

D15=010/(D2-Om011/(D5-09>  + 012  _ * _ __  

P16=YCSl37*01 

017=0CS13  7 

018=09*013 

P1?=D9»014 !_  

020=09*015 

021=016+018  _ _ _ 

D22=CS 1371 

023=021/(02-0170 +019/(D5-_D17)+D?0/(09iP17)+jD22 

024=020/(017-09) 

025=019/ ( D 17-05) 

' D 2 6= 02 1/(0 17 - 02 ) 

027=-02*T_  

D2o=-0  5*T 

0_29=-09^.r_ 

030=-017+f 

A H 7 1 37  =026  *EX  P(02 7)+0_2  5*  EX  P ( D 2 6 ) + 0 24 + EXP  ( D 29 ) + 0 23  + EXP(C3Q) 

ACT134=A9T 134*0CS134 

A CT 1 3 7=  A'l T 137+0CS137  . 

R=ACT 1 34/ ACT  1 37 

P PINT*  I =1! , PH  I , “ _T=1*iT,*^ _AMT134f**,AMT134,” A K T 1 3 7="tA  HT1  37 

U235I  = U235I*E  XP ( -C2+T ) 

PRItlT*»*HJ2_35I=",_y235I 

S B 1 3 3 1= C 5 * E X P ( -C  2 + T C 6 * E X P ( - C 3 *T  ) 

T El 33I=G13*EXP(- 010*7) +C14+EXP (-C2*T)-C15*EXP(-C3*T) _ 

I133I=C23+EXP (-o2*T ) -C24+ EXP ( - C3*T ) +C2 5*EXP (-C1 0 *T ) + 

1C26*EXP(-C20  + T)  _ 

XE133I =C3  5*"  EXP (-C2*T ) -C36+EXP (-C3*T) +C3 7 *EXP ( -Cl Q+T ) + 

103o*EXP(-C2Q*T) +C39+EXP (-C29+T)  

CS133I=C47+EXP(-C2*T) -C4o+EXP (-C3*T)+C49*EXP (-C10+T) + 

1C50*EXP(-C20*T) +C51  + EXP  (-C29M) +C52-EXP  (-C40  + T)_ 

CS134I =A9T 134 

I137.T=D7  + EXP( -02* T ) +D8*EXP(-D5+T)  _ 

XE 1 371=01 3*E X P (-02+T ) +0 1 4* EXP ( -05* T ) +015*EXP (-09 *T) 

CS 13  71  = A •IT  137  _ _ __  

101  CONTINUE 

STOP"ENO  OF  PROGRAM^  _ _ 

END  
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-n 

vO 

CD 

'JD 

rH 

-r 

ro 

ro 

ro 

m 

M 

• 

II 

* 

n 

• 

II 

• 

ii 

• 

ii 

• 

II 

• 

ii 

• 

li 

• 

li 

tM  PH 

C\J 

PH 

ro 

ph 

f\j 

PH 

PJ 

r-i 

P0 

ph 

ro 

PH 

ro 

r-H 

ro 

PH 

ro 

»-H 

ii 

ip 

ii 

in 

ii 

IT. 

n 

LP 

II 

in 

II 

in 

it 

in 

II 

in 

it 

in 

II 

H 

M 

ph 

ro 

PH 

ro 

PH 

ro 

PH 

ro 

*-H 

ro 

r-r 

ro 

»-H 

ro 

PH 

ro 

ro 

X 

<M 

X 

<vi 

X 

<vj 

r 

nj 

X 

cvj 

X 

CVJ 

X 

CVJ 

X 

CVJ 

X 

CVJ 

X 

CVJ 

CL 

3 

Q- 

3. 

a. 

r> 

a 

O 

CL 

z> 

a. 

ID 

a. 

ID 

a. 

ID 

CL 

ID 

CL 

i y 
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Program  YCSFOUR 


PROGRAM  YCSFOUR ( INPU T , OUTPUT) 

PEAL  11331,11371 

DIMEKSION  I KTFGR  ( o ) , 4R=>AY1<10)  , SPRAY  2 (1C)  , ARRAY3<4,2) 
DATA  INTEGR/3* (-0) ,0,2“ (-0)/ 

CALL  SYSrEMC(li!,INTEf,P) 

REA0M,U235I»SIGF,SIGA 
READ* ,SB1 33I,YS3133,0Sn133 
READ* ,TE1 33I,YTE133,3Trl33 
PEAO*'  ,H33I,YI133,jI13^ 

READ* , XR1 331, XE 1333, OXr 133 

READ* , CS1 331 , CS 1 33C 

READ* , II 37 1, Y 1 137, 01 137 

READ'  > YEl 371, DXE137 

PE AD* ,CS137T, YCS137,nCGl 37,CS13:3 

PF  AD'  , ARRAY1,  A RRAY2  , A RRA  Y3 

PRINT*  , ARRAY1 

PRINT*  , ARRAY? 

PRINT*, ARRAY3 
00  101  1=1,10 
PHI  = A RRAY 1(1) 

CO  101  J=l,  10 
T=ARPAY2v J) 

DO  101  K= 1, 2 
CS13<-  I=A?’RAY3(  1,*> 

YCS134=APRAY3(2,<> 

nCS134=ARRAY3(3,*) 

CS  13<*C=ARRA  Y3  (4,0 

C1A=U?35I*SI3F*PRI 

Cl  = YSBl 33*  Cl A 

C2=SI GA*°HI 

C3=DSB133 

C4=S?  1 331 

C5=C1 / (C3-C  2) 

C6= (C 1 ♦ <C2*  C4 ) - < 3 3* C 4 ) ) / (C3-C2) 

C7=YTE1 33* CIA 

Cfl=C3*C5 

C9=C3*C6 

C10=DTE133 

Cll=r7*CS 

C12=T  F133I 

ri3=Cll/(C?-C10)  ♦ Cl?-Cr)/(C7-ClO) 

C14=C11/(C10-C?) 
cis=rq/  (Cio-03) 

C16=YI131*C1A 
C17=f 1 0*C 1 3 
018=010*31^ 

C19=C10*C1G 

C?0=ril3? 

C21=C 16+C10 
C22=I 1 331 
C23=C  21/  ( C 20-C  2) 

C24=ri9/(C20-3.3) 

C2S=r 1 7/ (C  2 0-C  10) 

C?f>=C22*Gl7/<3  10-C20)  -0 19/ <C 3-C20) 

C27=yE133C*PHI 
C28=PXE13  3 
C29=f  27 +C23 
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74/74  OPT=  1 


FT N ‘».6*-h«,6 


( 


030=070*023 
C31=r20*024 
072=020*025 
C33=C  2 0*026 
C34=YF 1331 
C36=C  30/(C29-C?) 

036=0 31/(0 29 -03) 

037=032/(029-0  10) 

039  = 0 33/ ( C 29-020 ) 

C39=C  30/ (C2-C29)  -C 31  / ( 0 3-C 29)  *072/  ( 0 1 0 -C2 9)  «-C33/  (C20-C29)  *034 

C40=CS133C,‘PHI 

04  1=0  26*035 

042=028*036 

043=C  28*037 

044=0.26*033 

045=028*039 

C46=CS133I 

047=041/(0*0-0  2) 

04  8=0  4 2/  (C4  0-C  3) 

C49-C43/(C4C-C10) 

C50=044/(C40-C20> 

051=045/(040-0  29) 

C52=C41/(C2-C4  0)  -C 42/  ( C 3-C4  0 ) *043/  ( C 1 0 -C4  0)  *044/(020-040)  * 
1045/(029-0*0)  *046 
C53=YCS134*C1A 
C54=r si33C»P4I 
C55=CS134C,‘  DHI 
C56=OCS134 
057=054*047 


058=054*043 

059=054*049 

06  0=C54*C5  0 

C3i=r  54*C3  t 

062=054*052 

C53=C  S134 1 

064=055  *C5G 

065=053*057 

066=0  65/ ( C64-C  2) 

C67=-C2*T 

068=056/(054-03) 

CG9=-C3*T 

070=0  56/(064-0  10) 

C71=~C 1 0*  T 

C7  2=CC  0/ ( C64-C  20) 


C73=-C20*T 


074=061/(064-029) 

C75=-C29*  T 
076=062/(064-040) 

C77=-C40*T 

078=065/(02-064) -05 V (03-064) *059/ (01 0-t 64) *060/(020-064)* 
1061/(029-064)  *05  2/ (04  0 -064) *C63 
C?9=-C64*T 


AMT  134 =066* EXP  (0:7)  -C6ll*EXP(  C69>  *9  70*  EXP  (071)  *C72*EXP  (C73)  ♦ 
1074*FXP(C75)*C76,FXP(C77)*C76*EX3(C79) 
n*oiA 

02=02  , 8 ' 

P3=YI 137 


PROGRAM  YP^FOUR 


74/7* 


OPT  = l 


FT N 4.6*4**6 


04=01*03 
* \ 05=01137 

06=11371 
07=0**/ (OS  -02) 

oa=oe  *04/(02-05) 

O9=0XE137 

oio=P5’,n/ 

0ll=P5*0s 
D12=X  E 1 37 1 
ci3=nn/  (09-02) 

014=^11/ (09-n5) 

015  = 010/(02 -09)*D11/(0F- 09  )«-pi2 

016=YCS137*D1 

017=003137 

018=09-013 

019=f  9*014 
020=09*015 
021=016+013 
022=CS137I 

023=021/(02-017) *019/ (05-017) *070/ (09-017) *022 
024=02 0/(O17-D9) 

025=0 19/ ( D17-D5) 

020=021/(017-0?) 

027  = -02*  F 
023=-05*  T 
D29=-09*T 
- 03 D=-017V  T 

AMT 137=026*  EXP (027) + D2?  *EXF(028) +D?4+EXP( 029) +023*EXP(D30) 
ACT134=A9T13**mDC5l  34 
ACT 137  =AdT 1 37* OC 51 37 
R=ACT134/ACT137 

POINT  *f"PHI  T=",T,"  YC31  34="» YCS1 3~»"  AMT134=”,  AMT  134 

1"  R=",R 

191  CONTINUE 

ST OP**ENO  OF  PROGRAM" 

ENO 


O 
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' f ■ 


i-r 


. ■ 


- 


■***> 


l.E»13  2.E*13  7. <^17  i.cH3  *.E*13  b.E»l?  7.F*li  t.e*13  C.E*13  l.SM*. 

2592000.  513*0011.  77V  2 0 0 . 10"»7r000.  129i*'001.  1SF51000.  mbOOOr.  ’O/nnOO.  27'” 
0.  3.  1.07F-S  9 . . 000  003?  1.07-:-*  1.07--22 


PHl*l.£»t 3 
PHI*l.E»l 3 
®HI * 1 • E ♦ 1 7 
PHI * 1 • F ♦ 1 3 
PHI^1.E»1 3 
PHI *1 , E*  1 3 
PHI*1.E*1 3 
PHI 3 1 * F ♦ 1 ’ 
PHI*1.E*1 3 
ohI*1.E»i 3 
PHI*1,E*13 
PHI*l.F»li 
PHI*1.F»13 
PMI*1.EM  3 
PHI«l.E»l3 
PHJ ■ l « E ♦ 1 3 
PHI « 1 . F ♦ 1 3 
PHl»t.E*l 3 

PHI*l.p*l 3 
PhI*1.E»*.  3 
PHI*2.FH  3 
PHI*2.F*13 
PHI *2. E ♦ t 2 
PHI«?.e*H 
PHI*2.F»1» 
“HI*2.FH7 
PHI  * 2. E ♦! 3 


T»?i.:P?910  . 
T*I10’«019. 
T3?n  -Oil. 
T«77'r II". 
T377'’jni0. 
TslO^-OOl*'. 

T3117701C''. 

T3i:=5oir'. 
T3j  o ) in. 

Tsi'-; o ion . 
T = • Ve-Or  10*1 . 
T s * 4 1 I Q *>  1 0 . 
Til  3 1-  0 no. 
T i 20’  .0  m . 

',3?o7i.i  no. 
Ti  ?33  3 non. 

T3?’7, 31or , 

T=?rl’CTlfl. 

Ti?rP2on'i. 

'•7M30H, 

T3?om  n . 

Tsitoi  on. 
TiiiSi-HOO. 
T3»7»F0O'l. 
Ti 777?  **n. 
T*i')77nno. 


YO^tliiC.  5 MT  1 %*  1 7*.  • 9 c 57 Of  ?<*1  r i . 9 ?1 679 3 5£6 3" 1 
YC1131-  3.opni52  f ‘*T1343ii.i.i  362323073  Pi.  0 22SPOi  39  213*. T 
YTSlIiir.  4MT1  3U  = 55l».7Y  1*<»2213-  9 i . 0 5*.G  *77**559  5 1 

YFS13-  3.1009  05  2 17713*  *b7<) . 5 T5793  9 9 5 2 Fs  . C 55  243  95  *43  3 C~ 

YCS 1 3*-  * 0 . 14713731*17 .756372261  9 = . 0 <57 21  * 35C1.3 ?.  u 

YCE13* s.00010>2  !MTl-»i4S-lb?n.*.72bJai07  F 1 . 0 9549  if  5 2 2 7 7 4 > 

YC9i:<*,i3.  AMT13.37U3 .15*7277?-  9*  . 1 ?C  41  3549 15  ?n 

YFM7*.s.i)09JCi2  1 -r  t 7-3  29»-6  .*452?*, 274  I*  . 1 215  } 0917?  1 - > 
VC7 1 ?»• 3 o . H7T131  *-'>0*l.Efl,7??i-r  »*  . l*570-%*?7*52> 

YCS l’** .009395?  1 *T1  »-3<63E.i7«--i?99*  ■>  « . 1?  ’♦?  9 0 2“  Ob  3 « 

VCS 1 .2*.  = 0.  lit  13!  *<»i»9 .59l°2e9,i*  9*  . 131205 9375»77 

YF31J1.3.  00  1 105  2 »H’,l3’13fcoY0.20fcl  7>73u  2 » . 156  3*  C 3*7  3 

Y**SlTi. *0.  17T17:  >5  357  .o*.  941 95  37  9«  . 21t  5 »?  1 255 1 ? S 

YC5l»*3. 010005  2 l-T17-390’i*  .00151  3371.  93 . 21?3’7  7 253  3 : • 

YC3l.3-.30.  A7T17’  3HS6’  ,73«9*«?3  9«  . 2*7 1 ?C 4955 7 V 

YCSl’i.s.  0900  15  2 lHTt7'.«117»0. *0191777  9 « . 209  2 2 727  3 1 3 J 1 

YCSl’i.iO.  .17  7 1 7'-  1 1 ',6 3*  • S J£i* 0377  «3 . 27? ? 4697 1 0 3 57 

YCSl7‘is.O0O’J05?  iHTl3l.3l4F97. 97325712  9*  . >7 9 j 221 7*- ’ 

YCS13ii=0.  A7T  13-il7787.11**  7353  9=  . 30t  iS -o7  2 ’9  i 2 

YCSi. 2 1.*. 0931052  »-iT171.»17S<.5.15JU7J?3  » 1 . 3 099*  Cm  3 0 1 5 > 

YCS  1 7*  3 j , jHl  1 3-3577. 2^477  *6  588  93 . 0 4 3?1 3 7?39 59 / 

YCS13*  s. 0001952  tHTt  7t3Si,6.i7i.ViS30S2  F • . 04* > *0 21152 j? 2 
YT«17.  ifl.  nil  »!.3?i>V.7fc$70(,ie9  9i,  103165 5 957295 

Yr31  3*  = . OC  09  35  2 F 77  l 74 32*54 . 74,5*7  7755  f * . 1 95  3b  » 7 2o9fi2c 
YCS  1 3*  *0.  UTl  3*3«  2»».li.bi.a?633  9 ■ . 1 7*t  99  i 155  7f,  7 

YCS13*  *.n"019‘.»  SYTl74«f  3l4.47»rS3l>7  Ft . 175 377 7 32 7b 2* 
X*lS17‘-«0.  A7T17U311  770. 25007517  9«  . 24?  i 67 1 22  7 3 C 8 
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J&6-. 


* 


phi*3.e*i7  T«ie*"*i*ir>.  ycei’a*. 0000013  4*71*4*11435, o"»5.7}9<i2  9*.3*?ooj*6*376s 

*n  T*t?*?  p ii"  .#  vcsi’^o.  aitu^i^o*. 7043*373  >.= . 39caiJ53n*3* 

PMI*3.<>17  T»l®3r.OP91.  YESl*4  = . oiluOr- ’ AlTl  1*.=  ir?S9.r-  0=  " 3 ® **  2*  . 707la3668**t5 

PHI*?.F*!T  T*li55«101.  YC?1.7i*d.  ArtT17:.*2‘;47®.8o9457  0/  < r . 37  8 lit"  3®54  4 9 J 

PHI*2.E*11  T*i*,r»n.in<\  *:si*A=, 001005 ? 4*71*4*25567. *'C)75et>  9 * . '7  1?<  7 368  0 0 0 1 

phI*2.E»13  T*l5t*1  111.  y;Si*4*g.  4*713'  *34795. *08  23021  9*  . 4 3*35 1 93718t 

®HI*3.*>17  7*18.1.49101.  Y^Sl?**.  00  0 107  2 t*»13<  *34394.  t'Hll’i  s*.4Y*Il* 0694 773 

PHI*2.E*t*  r»7"7.M11.  YCtlTA  = 0.  4PTl*'-*<,*2* ..*o9s6£3  9 * . ‘ 95*.  3 4 3 17  91*4 

PHI*3.£*17  T.’>i!-lpv1p,  yesi***.000C05?  4 “T  1 3-  **4  73 1 . t-r  5 s 1 1 1 1 9*  . *5  6*9  36950  379 

PHI*3.E*15  ’=?*»*011".  YCS 1 74  * o.  AIT  13.  sEl  1 0 7 . 64  33o  H- 1 8*  . 561  24  544 3«  1 36 

PHI*3.£*1»  T****’!’!".  YCSlT4s.G0100i  3 A,'Tn-s5'?0o.7-7!i3‘.*  9 * . 56  7 3 3 05S*0  *91 

PHI*3.E*13  Tr?<-i->10in.  YCSl*<.*1.  4>1T  l 3- *6  5S6E  . 9k47«.39*  ’*  .61*  46a  pTd""?  Co 

"HI*2.e*l3  Ts??’'»1110.  YCSl.74*.  00011??  AMrt3A*60®81. 3*51791.  3 »«  . S 1 ▼ 0?  5 4-99  3 C 7 

PHI*  3. E*  t 7 Ts3rv*"in.  YCS!*‘=C.  A M7  17**1  1 1 ? . f.  ®E60859?  »■  . t 55019?  89*718 

PMI»3.C*17  T«353?"1P.  YC'IJi  *.  0000053  i M T 1 *4  = i 21 4 . g*  (,‘’47  7 9 7 K * . Out  3 *4  7 0 379 1 20 

PHl*3.E*l7  ’»5i«'0P1.  Y^rno.  1 -f  1 *.*5  4 7 4 . 1 5 3 7*  1 737  9=  . lb  374  55  1,0 1 f 5 

PHJ*3.EM  7 T = ?15i110.  YE313* =.1000032  4*T174=5a7?. 069)5*058  E* . 1 5*54 1 20  0797  3 

PHI=3.E»11  T*7.”6ni.  YCS1*'=C.  4*ri74*137*l. 74913344  9* . 3 6 31  3 94 id  1 75 3 

PHI*3.£»17  T*7’7? " 1 P . YCS17I  *.0001053  £871*4=13903.73717558  E * . 36  733  3 1 77  914  9 

PHI=3.E*l3  TalPVPVin.  YES 1 34  * 0 • «MT  1 7- a ’ .9  3 " . 37  34  1 1 39  7=  . 76 1 34298 109 3 3 

PHI  s 7.  E ♦ l 7 '.1  o*’111" . YCS174*.  0000053  4 “T  1 74*  350  0 1 . 1 7*  * 1 3 5 f *.  7633  it  3 354b  77 

PMI*3.F»13  T=1 355  P 1ll.  YCS134  = o.  4-1717.  = *17=i*. 33*3290.  9'  . 'tit  T 506755 1 1 3 

PHI-l.c  + lI  T=179<"010?.  Yf  Sl*4  = . POP  505  3 4 *T  1 74  = 38  53 .7. 9S79®  * 1 3 9*  . *4g3  0 37  99  77  17 

PHI  *3. £♦ 1 3 T=ly-;=10PP.  YCS1’4=0.  AH’  1 7-  = 5 ~ 19'  . 7455  01 3 ' 9* . 5541  73  9 154065 

PHI*  3. E ♦ 1 3 T=t5E-11V>.  yCS1?4  = . 00  0 C05  3 4 MT 1 3**55204 . 9*52 764 1 9=  . i>5?74  7 1 29940  3 

PHI*3.E»17  T= 1(141199.  *CS1*'.*G.  4 ‘i  1174*7757"". 4913117’  9* . 64" 7 7 8 4° 255 36 

PHI  * 3 . c ♦ 1 7 T*18tl  0 111.  YC'l  . 110005  3 4 MT 1 *4  *7  381 1 . 5 4 >*  74 1 = * . * 498  ^3  0 19  0741 

PHI*3.  = »15  T*  7074 1 1lf . YCS1?4*3.  tp T 1 71  * 3*35 » , t « 13 jg 9* . 7 41 759 Oo 9 39  — 

PHI*  3 . E*  1 3 T.70’41113.  ’rsi-".*. 900005  3 APT  1 3'*s94.  9 0 . 35?i8?5  1 9 = . 74345  •'lt.53549 

PHI  *3.  E*  1 3 T*  3777"110.  YCSlUsO.  tlT  1 7 . * 1 1 57  4 1 . 375;.7  ?•*  9*  . 3 71 748  4?94  3 77 

PHI  *3 . E * l * W:7’0JV'.  YC«"!?4*.  000  )053  4 - T 1 ya  * 1 If  i 39.  ’ 9 ' 792  d 9=  . 83290  3170371a 

PHI*3.E*13  T*»!  a?,'i9P.  ycsi*.=j.  AmT  17'-«  1407  (=8 .5P5  31.7'-  9*  . 92C  4575119  as  S 

PHI* 3.  E*  1 3 T* 3f  9 ' 1 )19.  yc' 1 3 4* . j 9 ) C05  2 4 -IT  l 1*08 co . °l  T595 4 9 *. 9 21492  330 75  Ja 

PHI**.  '"♦17  T=2:A',P10.  YES  1 * u . «H’i34s?n7.CE23t435i  -*.  08675757777153 

PHI*4.EM7  T*3‘  5*0  n."  YCE17- *.OC010C3  1 •*  Tt  74*  3 1 37 . *4’9C  1 2 3 S (■*.0079717***1476 

PHI*4.E*17  T*?l**t!-17.  YC517‘*0.  A -T  1 74  * 1 0 74?  ,4C  336f  5 1 9 * . 3 16*  3 S3  ^8  447  S 

PHI  *4  • E ♦ 1 7 T*5  10‘.  flPI.  YC513*  * .3000C1  3 f IT  1 34  *1  02*8  .7  1’5»Y5  b F *.  2 17775  97  339  0 3 

PHI*4.E*17  TsT’’*  941  • YC317(*0.  * H7  1 34  * 3-,  l 34  ,7f,l«  1 55 1 9 * . 749E  3 3a  1046  2e 

PH1*4.E*17  T=7'7',1on.  YE  S 1 7i  s . "0  9 0 Oa  ? 187  174*34217.4373389*  F* . 35  07»  0 9 3961 97 

PHI  ** , E M 3 T*197Y01,10.  YCS 174  = 0 • C * T 1 3- * 4 71 6 C . loi  5 Ob91  9*  . 48 1 54  » 3 07  3 0 0 2 

PHI  *4 . E » 1 T T*10*T()i(5«.  YES17L*. O0000S3  A *T  1 *4*  * * ?6  7.  ?P2  3 73  7 1 9=  . 4 3 35a  94 16338  2 

PHI  * * . E ♦ 1 3 7*13*59  Oil.  Y<;si3**0.  AlTl* -*5657*. 7 98  754  9.  . j 1 121  32  7 9o097 

PHI  *4 . * ♦ i 7 7*  1390071 0.  Y2S:?4  = .90  J0052  4 -T 1 3a*  Oc.76  9 . 4257  7?  £ , 9 * . 6 1 273  7 39^  8 17  7 

°HI *4  « E ♦ l 7 T*lr""0  711.  ’OElJtrO.  4*T17 -*C*14’.0295  297Y  2* . 7 78 45 354773 44 

®HI*4.  E*  1 3 Til'iO-lllo.  YCS134*.  00  109*  3 48713**  94  38  3. £5 15  281*  9*  . 7 39*o  e 8 9 5*9  71 

PHI *4 . E ’ l 7 Ti’I’^OP^P.  YC313**P.  4<T13' *127172. r»6711Y  9* . 8631355177508 

PHI  *4  . E » 1 3 T*l.«l<- 0 70".  YCE13=  *.  OOOlOa’  4 1’  1 3*  * 12c  3 09 . 155  91 3 3 9 * . »64»  1 93  3 38  1 34 

PHI*4.E*13  T*20’-1901.  YES134*1.  AiTl*'  *1  •190.'*  0.301  71 1 6 9*  . ‘*3*  51  j 7 8 7 1 5 ? 1 

PHI *4 • E ♦ t 3 T*?9’4  3 IIP.  YES  1 34  * . 00 0 00*  7 4 HT 1 3« * 159241 .25507 8*  2* . 56o 5 7*  1 9436® 7 

PHI  *4  . E * 1 3 r*37*70110.  v-'Sl?4*o.  *.PTt7-*t  15469.478*041  2*  t . 10  49  0 0T425  48 

PHI  *4.  E*  1 7 T*?”7m".  YC51  74*. 000115  3 4*T  l 3"-*  l9?o3  3 . 18®  593  0 9*1.1359  38  7 07  545 

PHI**.  **t7  r*?*9?9P9J.  yeeiioso.  A*T  1 34  *2  74  0 46  . gi.g  67  ?’  **1.2314  376*5996 
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